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ABSTRACT 


The  minimum  resolvable  temperature  difference  (MRTD)  is 
widely  accepted  as  the  parameter  that  best  describes  the  field 
performance  of  a  thermal  imaging  system  (TIS) ,  Mathematical 
modeling  that  accurately  predicts  the  MRTD  has  been  of  major 
interest  to  the  infrared  community  over  the  last  30  years. 
This  work  reviews  the  currently  accepted  models  for  predicting 
the  MRTD.  Simplifying  assumptions  used  by  these  models  which 
deal  with  the  target  spectrum  are  discussed  and  tested  using 
specifications  taken  from  a  standard  forward  looking  infrared 
(FLIR)  system.  In  addition  new  models  are  proposed  and 
tested.  Two  of  these  models  are  a  direct  extension  of  the 
recently  proposed  Vortman-Bar-Lev  adaptive  matched  filter.  A 
third  model  is  based  on  the  novel  concept  that  the  MRTD  curve 
is  predictable  from  a  threshold  condition  on  the  visibility, 
rather  than  the  signal-to-noise  ratio,  of  the  system-degraded 
bar  pattern. 


•  •  • 

111 


Accession  Per 

NTIS  GRA&I 
DTIC  TAB 
Unannounced 
Justification. 


I  By - 

Jibuti  on/ 


Availability  Codes 
(Avail  and/or 
Wet  j  Special 


□  □ 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  .  1 

A.  REVIEW  OF  LITERATURE  .  1 

B.  OVERVIEW  OF  THE  THESIS .  2 

C.  FUNDAMENTALS  OF  THERMAL  IMAGING  SYSTEMS  ...  3 

D.  BASIC  PARAMETERS  .  5 

E.  MINIMUM  RESOLVABLE  TEMPERATURE  DIFFERENCE 

(MRTD)  .  9 

F.  NECESSITY  OF  A  MODEL .  10 

G.  MOTIVATION  FOR  THIS  WORK .  12 

II.  ESTABLISHED  MRTD  MODELS .  16 

A.  SYSTEM  DESCRIPTION  .  16 


1.  NETD  Calculation  for  the  Sample  System  ...  17 

B.  THE  RATCHES  MODEL  FOR  THE  SAMPLE  SYSTEM  ....  20 

C.  THE  LLOYD  APPROXIMATE  MODEL  FOR  THE  SAMPLE 

SYSTEM .  22 

D.  THE  VORTMAN  AND  BAR-LEV  MODEL  FOR  THE  SAMPLE  SYS¬ 
TEM  .  25 

III.  PROPOSED  MRTD  MODELS .  30 

A.  MODIFIED  VORTMAN  AND  BAR-LEV  MODEL  FOR  THE  SAMPLE 

SYSTEM .  30 


iv 


B.  A  VISIBILITY  MODEL  FOR  MRTD  MEASUREMENTS  ...  33 

C.  INCLUSION  OF  SAMPLING  ARTIFACTS.  EFFECTS  ON  THE 

RATCHES  MODEL  .  41 

IV.  COMPARATIVE  ANALYSIS  AND  CONCLUSIONS  .  47 

APPENDIX  A  -  INTERPRETATION  OF  THE  RATCHES  MODEL  .  .  50 

A.  PRELIMINARIES .  51 

B.  DERIVATION  OF  THE  NOISE  EQUIVALENT  TEMPERATURE 

DIFFERENCE  (NETD)  .  54 

C.  DERIVATION  OF  THE  MINIMUM  RESOLVABLE  TEMPERATURE 

DIFFERENCE  (MRTD)  .  65 

APPENDIX  B  -  DESCRIPTIVE  FLOWCHARTS  OF  SIMULATION  ...  79 

APPENDIX  C  -  EXACT  FOUR-BAR  FOURIER  SPECTRUM  .  80 

APPENDIX  D  -  SAMPLE  SYSTEM .  87 

LIST  OF  REFERENCES .  89 

INITIAL  DISTRIBUTION  LIST  .  93 


V 


TABLE  OF  SYMBOLS 


a:  Modified  visibility  parameter 

/9:  Parameter  to  express  ATSC  as  a  function  of  NETD 
T:  Constant  to  convert  units  of  noise  [°C2  mrad] 

A 0(A):  Wavelength  dependent  part  of  P(A  t)  [watts/^m] 

Afn:  Noise  equivalent  electrical  bandwidth  [Hz] 

Af^:  Noise  equivalent  electrical  bandwidth  [cy/mrad] 

AT:  Target-background  temperature  difference  [°C] 

ATS:  Apparent  temperature  difference  (Mod.  Visib.)  [*C] 
AT8C:  Threshold  ATS  for  resolution  (Mod.  Visib.)  [°C] 

Ax:  Horizontal  DAS  (mrad] 

Ay:  Vertical  DAS  [mrad] 

Ayf :  Angular  distance  between  scan  lines  (Ay/r?ovsc)  [mrad] 
r?0(A):  Optical  efficiency  of  the  viewer 
rj wsc;  Overscan  Ratio  for  the  device 
tjsc:  Scan  efficiency 

A.:  Wavelength  of  monochromatic  infrared  radiation  [Mm] 
Ap:  Characteristic  wavelength  of  the  detectors  [Mm] 
a:  Root  mean  square  noise  at  the  output  of  [volts] 

7:  Threshold  signal-to-noise  ratio  for  resolution 
a:  CRT  spot  size  parameter 

A:  Normalized  spatial  frequency  for  optics  MTF  [Hz] 
Ad:  Area  of  a  single  detector  [cm2] 

vi 


At:  Total  Area  of  the  target  [mrad2] 
bj(x):  Noise  along  the  ith  video  line 

d:  Focal  length  of  collecting  lens  [cm] 

D:  Diameter  of  collecting  lens  [cm] 

D*(A):  Specific  detectivity  as  a  function  of  wavelength 
[cm (Hz)  -5/watt] 

DAS:  Detector  Angular  Subtense  [mrad] 

F:  F/  Number=d/D 
FLIR:  Forward  Looking  InfraRed 

fox:  Horizontal  spatial  frequency  of  the  target  [cy/mrad] 
fr:  3-dB  frequency  of  electronic  Roll-Off  [Hz] 

Fr:  Frame  rate  [Hz] 

frx:  3-dB  frequency  of  electronic  Roll-Off  [cy/mrad] 
fx:  Spatial  frequency  along  the  horizontal  direction 
[cy/mrad] 

fy:  Spatial  frequency  along  the  vertical  direction 
[cy/mrad] 

H(a):  Spectral  irradiance  onto  a  detector  [watt/ (cm2  /xm)  ] 
Hd:  Display  MTF 

Hi,:  Overall  System  MTF  ( HoptHdetHELECTH(J) 

Hdet:  Detector  MTF 

hoET-ELECT:  Detector-electronics  impulse  response  [volts/watt] 
haECT:  Electronics  impulse  response 
HElect;  Electronics  MTF 
HFOV:  Horizontal  Field  of  View  [mrad] 

Hmf:  Matched  Filter  MTF 


vii 


V 


^"(x,y)  * 
^T(x,y)  * 


I 


T* 


k: 


K 

L 

JCA) 


M 


<*,y) 
MRTD 

MTF 

n(x.y) 

NEPu> 

NETD 

NETD' 

np 

PU,t> 


PSD: 


r: 


R: 


Optics  MTF 

Normalized  Fourier  transform  of  the  target 
Normalized  Fourier  transform  of  one  bar  (along  x) 
Normalized  time  dependent  part  of  P(A  t) 

1{t)*r(t)^ELECT<t) 

Spatial  description  of  the  displayed  target 
Spatial  description  of  the  actual  target 
Fourier  Transform  of  the  target 
constant  to  transform  temperature  into  radiant 
displayed  energy  [watt/*C] 

Constant  to  define  a  Matched  Filter 
Length  of  the  bars  of  the  standard  target  [mrad] 
Spectral  Radiance  from  source  (watts/ (cm2  st  /xm)  ] 
Signal  energy  per  unit  area  of  display  [  Joule/mrad2] 
Minimum  Resolvable  Temperature  Difference  [*C] 
Modulation  Transfer  Function  (all  are  normalized) 
Spatial  description  of  displayed  noise  [Joule/mrad2] 
Spectral  Noise  Equivalent  Power  [watts] 

Noise  Equivalent  Temperature  Difference  [*C] 

NETD  converted  into  radiant  displayed  energy  [Joule] 
Number  of  detectors  in  parallel 

Time-varying  spectral  radiant  flux  on  a  detector 
[watts/ Mm] 

Power  Spectral  Density  [Watts/Hz] 

Detector  response  function  [volts/watt] 
Autocorrelation  function 


viii 


R(a):  Wavelength  dependent  part  of  the  Fourier  transform 
of  r 


R(f):  Frequency  dependent  part  of  the  Fourier  transform  of 
r 

tt:  Eye  integration  time  [sec] 

S '  j :  Single-sided  Noise  PSD  at  system  input  [Watts/freq] 
S'0:  Single-sided  Noise  PSD  at  system  output  [Watts/freq] 
Ss:  Double-sided  Noise  PSD  at  system  input  [Watts/freq] 
S0:  Double-sided  Noise  PSD  at  system  output  [Watts/freq] 
t:  time  [sec] 

T:  Temperature  [*C] 

Tb:  Background  Temperature  [ * C] 

TIS:  Thermal  Imaging  System 
VFOV:  Vertical  Field  of  View  [mrad] 

Vn:  Noise  Voltage  at  the  output  of  HELECT 
V#:  Signal  Voltage  at  the  output  of  HELECT 
v#(t):  Instantaneous  output  voltage  from  detector-electron¬ 
ics  [volts] 

Td:  Detector  dwelltime  [sec] 

vx:  Horizontal  scanning  velocity  [mrad/sec] 

W:  width  of  a  single  bar  of  the  standard  target  [mrad] 
y. :  Vertical  coordinate  of  the  ith  video  line  [mrad] 


ix 


J 


ACKNOWLEDGEMENTS 


I  would  like  to  thank  Professor  Ron  J.  Pieper  for  his 
continuous  guidance  in  the  pursuit  of  this  thesis.  It  is  also 
my  sincere  duty  to  acknowledge  my  EW  Curricular  Officer  and 
his  staff  for  their  support. 

I  would  like  to  thank  the  Argentine  Navy  and  the  U.S.  Navy 
for  providing  this  opportunity. 

Lastly  I  would  like  to  express  my  sincere  thanks  to  my 
wife,  Martha,  for  her  support  and  patience  during  my  stay 
here. 


x 


I.  INTRODUCTION 


A.  REVIEW  OF  LITERATURE 

The  minimum  resolvable  temperature  difference  (MRTD)  has 
been  widely  accepted  as  the  parameter  that  best  describes  the 
field  performance  of  a  thermal  imaging  system  (TIS)  when 
involved  in  recognition  and  identification  tasks  [Ref.  1], 
[Ref.  2],  [Ref.  3],  [Ref.  4]. 

It  became  a  popular  design  tool  for  forward  looking  infr¬ 
ared  (FLIR)  systems  in  1975,  when  the  U.S.  Army  Night  Vision 
Laboratory  (NVL)  published  a  successful  mathematical  model 
[Ref.  5]  which  was  based  on  matched  filter  theory  as  an 
alternative  to  the  predated  perfect  synchronous  integrator 
model  [Ref.  6].  In  spite  of  its  success,  shortcomings  have 
been  encountered  in  both  mathematical  modelling,  and  measure¬ 
ment  of  the  MRTD  for  TIS  [Ref.  7],  [Ref.  8]. 

To  improve  the  MRTD  prediction,  other  models  have  been 
proposed,  by  either  modifying  the  original  from  Ratches 
[Ref.  8],  including  particular  concepts  for  specific  types  of 
systems,  or  using  generic  expressions  valid  for  any  type  of 
TIS  [Ref.  9]. 

Perhaps  the  most  critical  issues  in  MRTD  modelling  are: 

1.  visual  perception,  which  has  been  an  area  of  fundamental 
research  by  several  authors  such  as  Johnson,  Blackwell, 
Schade,  Genoud,  Sendall,  Rose  [Ref.  10],  Rosell,  Willson 
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[Ref.  11],  Kornfeld  and  Lawson  [Ref.  12],  [Ref. 13], 
[Ref.  14],  [Ref.  15],  [Ref. 16]. 

2.  effects  of  sampling,  [Ref.  17],  [Ref.  18],  [Ref.  19], 
[Ref.  20],  [Ref.  21]. 

3.  noise  characteristics/  [Ref.  5],  [Ref.  8],  [Ref.  2]. 

With  respect  to  MRTD  measurements,  objectivity  has  been 

the  main  goal.  Several  laboratories  are  currently  offering 
high  quality  MRTD  measurements  following  procedures  that, 
though  very  similar  in  principles,  slightly  vary  according  to 
special  needs.  The  use  of  computers  and  voltage  prior  to  the 
display  to  eliminate  the  logical  operators'  subjectivity  is 
being  experimented  using  not  only  vertical  targets  but  also 
rotated  ones  for  paarticular  applications  [Ref.  1],  [Ref.  2], 
[Ref.  23],  [Ref.  24],  [Ref.  25],  [Ref.  26]. 

B.  OVERVIEW  OF  THE  THESIS 

The  present  chapter  introduces  the  fundamentals  of  the 
MRTD  as  a  figure  of  merit  for  a  TIS  and  the  discrepancies 
currently  reported  between  predicted  and  measured  data. 

Chapter  II  briefly  describes  MRTD  models  from  Ratches 
[Ref.  5],  Lloyd  [Ref.  10,  p.  184],  and  Vortman  and  Bar-Lev 
[Ref.  8]  that  apply  matched  filter  theory  and  presents  their 
predictions  for  the  sample  system  given  in  Appendix  D. 

Chapter  III  considers  four  proposed  models:  the  first  two 
modify  the  matched  filter  presented  by  Vortman  and  Bar-Lev, 
the  third  uses  the  concept  of  visibility  of  the  system-degra- 
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ded  bar  pattern,  and  the  fourth  shows  the  effects  of  sampling 
artifacts  on  the  original  Ratches  model. 

In  Chapter  IV  a  comparative  analysis  between  the  different 
models  is  performed,  their  feasibility  analyzed,  and  the  final 
conclusions  stated. 

Appendix  A  covers  the  development  of  Ratches  model , 
Appendix  B  presents  a  flowchart  to  provide  a  pedagogical 
description  of  the  simulations  performed  in  MATLAB,  Appendix 
C  includes  the  spectral  analysis  of  the  standard  MRTD  target, 
and  Appendix  D  lists  the  data  of  the  sample  system  used  for 
the  different  simulations. 

For  typographical  convenience  a  non-conventional  notation 
has  been  adopted  for  all  the  equations  in  this  work.  Specifi¬ 
cally,  the  functional  dependence  is  shown  as  a  subscripted 
argument  leading  to  a  more  compact  form,  i.e.,  h1(x,y)  is 
represented  as  hj(xy). 

C.  FUNDAMENTALS  OF  THERMAL  IMAGING  SYSTEMS 

An  object  that  presents  an  apparent  temperature  difference 
with  respect  to  the  background  can  be  detected,  and  subse¬ 
quently  recognized  and  identified  by  means  of  a  thermal 
imaging  system.  This  difference  in  apparent  temperature  has 
to  be  large  enough  to  overcome  variations  in  the  background 
temperature  and  other  sources  of  noise,  as  well  as  the  even¬ 
tual  attenuation  produced  by  the  atmosphere  along  the  energy 
path.  The  term  apparent  describes  the  fact  that  different 
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emissiviti.es  may  need  to  be  considered  for  target  and  back¬ 
ground  . 

Due  to  their  popularity  and  simplicity,  only  conventional 
scanning  thermal  imaging  systems,  known  as  FLIRs,  will  be 
considered  within  the  scope  of  this  work.  Consideration  of 
other  recent  designs,  e.g.,  mosaic  staring  infrared  systems, 
infrared  linescanners,  or  pyroelectric  array  systems  that 
present  particular  characteristics,  would  require  an  extension 
of  this  investigation  beyond  the  limitation  in  time  imposed. 
This  section  follows  the  description  presented  by  Lloyd 
[Ref.  10,  p.  15],  whose  text  has  been  accepted  as  a  fundamen¬ 
tal  reference  by  the  entire  infrared  community.  Figure  1.1 
shows  one  possible  way  of  implementing  a  FLIR. 

The  optical  system  collects,  spectrally  filters,  spa¬ 
tially  filters,  and  focusses  the  radiation  pattern  from 
the  scene  onto  a  focal  plane  containing  a  single  small 
element.  An  opto-mechanical  scanner  consisting  of  a  set 
of  two  scanning  mirrors,  one  sweeping  vertically  and  the 
other  horizontally,  is  interposed  between  the  optical 
system  and  the  detector.  The  ray  bundle  reaching  the 
detector  from  the  object  moves  as  the  mirrors  move, 
tracing  out  a  TV-like  raster  in  object  space  as  shown  in 
Figure  1.1.  This  process  of  detecting  the  scene  sequen¬ 
tially  is  called  scene  dissection  [Ref.  10,  p.  8]. 

The  energy  of  the  electromagnetic  field  incident  onto  the 
detector  produces  a  response  in  the  form  of  an  electronic 
signal  that  is  then  processed  by  amplification  and  filtering 
circuits. 

Finally,  the  TIS  requires  a  display  in  synchronism  with 
the  scanning  components,  such  that  the  image  is  presented  to 
the  observer  in  an  appropriate  way. 
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Detector  projection 


Vertical  Scanner 


Display 


Figure  1.1.  Simplified  FLIR  (After  Ref.  10,  p.  10) 

The  system  can  be  represented  by  a  block  diagram  of  four 
basic  cascaded  components,  which  filter  the  incoming  electro¬ 
magnetic  signal:  optics,  detector,  electronics,  and  display. 
Figure  1.2  depicts  this  concept  and  introduces  the  symbology 
of  the  individual  transfer  functions. 

D.  BASIC  PARAMETERS 

A  system  with  such  a  diversity  of  components  as  a  TIS  has 
a  great  number  of  parameters  that  describe  the  particular 
performance  of  each  element.  Braddick  and  Ludlow  [Ref.  1] 
list  more  than  twenty  different  parameters  of  importance  in 
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Optics  Detector  Electronics  Display 


Figure  1.2.  Block  Diagram  of  Simplified  FLIR 

TIS  design  and  evaluation.  However,  finding  a  single  parame¬ 
ter  to  describe  the  overall  performance  of  a  system  has  always 
been  of  primary  interest  for  engineers. 

From  the  beginning,  the  infrared  community  identified  two 
basic  aspects  of  fundamental  importance  in  TIS  assessment: 
thermal  resolution  and  spatial  resolution.  Historically,  they 
lead  to  the  first  two  parameters  of  common  use:  the  Noise 
Equivalent  Temperature  Difference  (NETD)  and  the  Detector 
Angular  Subtense  (DAS) .  The  former  is  related  to  thermal 
resolution;  the  latter,  to  spatial  resolution. 

The  NETD  can  be  considered  as  "a  measure  of  the  ability  of 
a  system  to  discriminate  small  signals  in  noise"  although  a 
more  proper  definition  follows:  "the  NETD  is  the  blackbody 
target-to-background  temperature  difference  in  a  standard  test 
pattern  which  produces  a  peak-signal  to  rms-noise  ratio  (SNR) 
of  one  at  the  output  of  a  reference  electronic  filter  when  the 
system  views  the  test  pattern".  (Ref.  10,  p.  166) 
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If  the  temperature  difference  between  target  and  back¬ 
ground  is  denoted  by  AT  (the  target  must  be  at  a  higher 
temperature  than  the  background) ,  the  signal  voltage  by  Vs  and 
the  rms  noise  voltage  by  Vn,  the  NETD  can  be  obtained  from: 

AT.V- 

NETD  =  _ -  .  (1.1) 

V 

vs 

It  follows  that,  if  the  signal  voltage  equals  the  noise  volt¬ 
age,  the  NETD  is  the  temperature  difference  between  target  and 
background.  Figure  1.3  shows  the  general  appearance  of  the 
target-background  arrange  as  well  as  a  typical  voltage 
waveform  for  a  scanned  line.  The  observed  DC  droop  is  an 
unwanted  signal  degradation  due  to  the  time  constant  of  the 
dc-blocking  circuitry  employed  to  couple  the  detector  and  the 
required  amplifying  electronics  and  can  be  minimized 
[Ref.  10,  pp.  336-345]. 

Equation  (1.1)  is  a  basic  expression  that  describes  the 
concept  of  the  NETD  as  a  figure  of  merit.  Appendix  A  provides 
a  detailed  explanation  of  the  measurement  process  and  a 
mathematical  model  for  the  NETD. 

On  the  other  hand,  the  more  intuitive  concept  of  DAS  can 
be  defined  as  the  angle  subtended  by  the  projection  of  the 
limits  of  a  detector  on  the  object  plane.  It  is  customary  to 
denote  the  horizontal  DAS  by  Ax,  and  the  vertical  DAS  by  Ay. 
The  time  required  to  scan  through  an  horizontal  angular 
subtense  is  known  as  dwelltime.  Figure  1.4  shows  a  magnified 
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Figure  1.3.  NETD  target  and  voltage  waveform 
(After  Ref.  10,  p.  167) 

projection  of  the  DAS  on  the  object  plane,  as  well  as  the 
horizontal  (HFOV)  and  vertical  (VFOV)  fields  of  view. 

But  this  information  is  not  descriptive  enough.  Thermal 
Imagers  have  been  mainly  used  to  recognize  and  identify 
"targets".  Therefore,  it  is  reasonable  to  look  for  a  pa¬ 
rameter  that  describes  their  overall  performance  in  those 
terms.  As  it  is  discussed  in  Appendix  A,  the  NETD  and  DAS  are 
employed  in  the  determination  of  a  more  elaborate  and  appro¬ 
priate  parameter:  The  Minimum  Resolvable  Temperature  Differ¬ 
ence  (MRTD) . 
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Figure  1.4.  Detector  Angular  Subtense  (DAS)  (After  Ref.  10) 

E.  MINIMUM  RESOLVABLE  TEMPERATURE  DIFFERENCE  (MRTD) 

The  results  obtained  by  Schade  [Ref.  10,  p.  183]  for 
resolution  in  photographic,  motion  picture,  and  television 
systems  were  adapted  by  Genoud  and  Sendall  [Ref.  10,  p.  183] 
to  infrared  systems. 

Infrared  imaging  systems  are  commonly  involved  in  recogni¬ 
tion  and  identification  of  military  ground  vehicles.  It  was 
shown  [Ref.  27]  that  the  probability  of  recognizing  and 
identifying  a  main  battle  tank  is  related  to  the  ability  of 
the  system  to  resolve  bar  targets  in  a  laboratory  experiment. 


The  infrared  community  has  selected  the  MRTD  test  to  be  the 
industry  standard  for  characterizing  the  performance  of  a 
thermal  imaging  system.  The  standard  pattern  was  chosen  to  be 
consisting  of  four  vertical  bars  of  7:1  aspect  ratio  (the  bar 
height  is  seven  times  the  bar  width)  on  a  background  of 
uniform  temperature.  The  bar  centers  are  separated  by  twice 
the  width  of  the  bars.  Currently,  most  MRTD  measurements  are 
performed  using  this  target. 

The  difference  in  apparent  temperature  between  the  target 
and  the  background  is  controlled  until  the  pattern  becomes 
recognizable.  (If  human  observers  are  employed,  majority  of 
opinions  is  required  to  accept  that  the  target  was  recog¬ 
nized.)  That  temperature  difference  which  brings  the  signal 
above  the  noise  and  makes  the  four  bars  visible  is  the  value 
of  MRTD  for  the  particular  spatial  frequency  of  the  target. 
Figure  1.5  shows  a  sample  set  of  four  targets  of  different 
spatial  frequencies.  It  is  noted  that  as  the  spatial  frequen¬ 
cy  increases,  the  overall  target  size  decreases. 

F.  NECESSITY  OF  A  MODEL 

The  need  for  a  mathematical  model  to  represent  the  MRTD  is 
quite  obvious  since  it  is  the  best  tool  a  designer  can  count 
on  to  forecast  the  task  performance  of  a  future  infrared 
imaging  system.  Unfortunately,  in  the  process  of  recognition 
and  identification,  the  ability  of  the  observer's  eye  and 
brain  to  handle  the  image  presented  on  the  display  is  of 
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Figure  1.5.  Sample  four  different  spatial  frequency  targets 


fundamental  importance,  and  then  the  model  must  include  a 
complex  [Ref.  13]  and  not  completely  understood  [Ref.  7] 
process  involving  eye  and  brain. 

The  Ratches  model  [Ref.  5]  defines  the  MRTD  in  terms  of  a 
scanning  system's  elemental  parameters  and  represents  the 
human  eye-brain  system  by  a  matched  filter.  It  is  extensively 
treated  in  Appendix  A  and  is  the  basis  of  this  work.  A 
typical  predicted  MRTD  curve,  corresponding  to  the  sample 
system  (Appendix  D)  that  is  employed  in  this  investigation,  is 
shown  in  Figure  1.6  (approximated). 


Predicted  MRTD  (Sample  FLIR) 

3  14- 


0  .1  .2  .3  .4  .5  .8  .7  .8 

spatial  frequency  (Cy/nRad) 

Figure  1.6.  Sample  MRTD  curve 

G.  MOTIVATION  FOR  THIS  WORK 

In  spite  of  the  success  of  the  Ratches  model,  there  still 
is  discrepancy  between  predicted  and  measured  data.  Several 
authors  have  reported  underestimation  of  the  MRTD  in  the  low 
spatial  frequencies  (too  optimistic  prediction)  as  well  as 
differences  in  the  high  spatial  frequencies.  Figures  1.7  to 
I. 10  clearly  show  this  situation  as  presented  by  three  diffe¬ 
rent  authors  for  four  typical  systems.  Figures  1.7  and  1.8, 
from  Vortman  and  Bar-Lev  [Ref.  8),  include,  besides  the 
Ratches  model  (represented  by  RL)  and  measured  data,  two  other 
models:  LSI,  not  discussed  in  this  thesis,  and  AMF,  presented 
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in  Chapter  II.  Both  frequency  and  temperature  are  normalized 
with  respect  to  unknown  (to  this  author)  constants  fg  and  K,  2, 
respectively.  In  both  figures  a  general  optimistic  behavior 
of  the  Ratches  model  is  observed.  Figure  1.9  from  Braddick 
and  Ludlow  [Ref.  1]  shows  a  calculated  MRTD  based  on  a 
different,  though  similar,  model  [Ref.  24],  and  Figure  I. 10 
from  McCracken  and  Wajsfelner  [Ref.  7]  plots  predictions  from 
the  Ratches  model.  These  two  figures  show  a  mixed  behavior  of 
the  Ratches  model:  too  optimistic  in  the  low  spatial  frequen¬ 
cies  and  too  optimistic  in  the  high  spatial  frequencies  when 
compared  with  measured  data. 


Figure  1.7.  MRTD  vs.  Spatial  Frequency 
(from  Ref.  8) 
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Figure  1.8.  MRTD  vs.  Spatial  Frequency 
(from  Ref.  8) 
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Figure  1.9.  MRTD  vs.  Spatial  Frequency 
(from  Ref.  1) 
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II.  ESTABLISHED  MRTD  MODELS 


A.  SYSTEM  DESCRIPTION 

Only  the  effects  of  the  four  cascaded  transfer  functions 
presented  in  Figure  1.2  have  been  considered  to  obtain  the 
overall  transfer  function  of  the  system.  Expressed  in  hori¬ 
zontal  spatial  frequency  domain  (fx)  ,  these  are  [Ref.  28]: 

1.  Optical  diffraction  limit  (H^)  : 

H0pt(fx)  =  ~  { cos'1  (A)  -A)] (1-A2)  }  <II*1> 


where  A=(lFfx/d)  1000,  F  is  the  f /number  and  d  is  the 
focal  length  using  the  units  defined  in  Appendix  D.  This 
strictly  positive  function  is  nonzero  over  the  range  of 
the  normalized  spatial  frequency,  A,  equal  to  [0,1],  It 
is  derived  from  the  autocorrelation  of  the  pupil  function 
of  a  spherical  lens  for  an  incoherent  optical  system. 
[Ref.  10,  pp.  99-103] 


2.  Detector  spatial  filter  (H^)  : 

_sin(xfxAx) 
“  (7rfxAx) 


(II. 2) 


where  Ax  is  the  detector  angular  subtense  (x  direction) . 

3.  Electronics  (HELECT) : 


HELECT<fx)- 


N 


(II. 3) 


The  MTF  Helect  is  more  commonly  expressed  as  a  function  of 
temporal  frequency.  The  conversion  from  temporal  to 
spatial  frequency  can  be  performed  by  the  relation: 
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(II. 5) 


f 


x 


f 


V 


x 


where  fx  is  the  spatial  frequency,  f  is  the  temporal 
frequency  and  vx  is  the  horizontal  angular  scanning 
velocity  [mrad/sec] .  [Ref.  5],  [Ref.  9],  [Ref.  28] 


4.  CRT  display  or  monitor  (Hd) 


(-afx> 

Hd«x>=e 


(II. 4) 


This  is  equivalent  to  assuming  that  the  display  can  be 
described  by  an  incoherent  point  spread  function  with  a 
gaussian  space  dependence  [Ref.  29]. 

Figure  II. 1  shows  the  four  individual  modulation  transfer 
functions  and  Figure  II.  2  shows  the  overall  modulation 
transfer  function  of  the  system. 

1.  NETD  Calculation  for  the  Sample  System 

The  complete  derivation  of  the  NETD  and  MRTD  models 
presented  by  the  Night  Vision  Laboratory  are  included  in 
Appendix  A,  where  the  NETD  expression  is  repeated  here  for 
reader  convenience  (A. 32): 


NETD  = 


4F2(Afn)  2 


3L 


(A) 


D*  dA. 


~5T  M 


(II. 6) 


Several  standard  approximations  are  applied  to  (II. 6)  in  order 
to  facilitate  the  NETD  calculation.  By  linearizing  the 
specific  detectivity  (D*(x)  =  AD*(>p)/lp) ,  and  considering  the 
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Figure  II. 1.  Individual  MTFs  of  the  sample  system 
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Figure  II. 2:  Overall  MTF  of  the  sample  system 
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actual  spectral  bandpass  of  the  system  (8-11.5  /xm)  ,  the 
integral  in  (II. 6)  can  be  replaced  by: 


.  rll.5  5L<a)  >. 

^o(ap)%>J8  stt* 


(II. 7. a) 


Employing  the  approximation  from  Lloyd  [Ref.  10,  p.  174], 
based  on  the  radiation  slide  rule, 


— ^?-!tb“--2LU,tb>  (II. 7. b) 
A.Tb 

where  c2=l.  4388xl04  [nm  *C]  and  TB  is  the  background  tempera¬ 
ture.  The  evaluation  of  the  integral  in  (II. 7. a),  for  the 
sample  system  yields  4.71xl0'5  [watt/ (cm2  st*C)]  and  the  ap¬ 
proximate  NETD  (II. 6)  is  then  evaluated  using  specified  values 
taken  from  Appendix  D: 


NETD 

= _ 4(2 . 5)  V29400  [HZ]  _ 

ir[st]  (0.005[cm])  (0.8)  (2xlO,0[  cm...H£-..  ])  (4 . 71xl0'5[ ...  ]  ) 

watt  cm2st  *  C 


-0.36’C 


B.  THE  RATCHE8  MODEL  FOR  THE  SAMPLE  SYSTEM 

As  derived  in  Appendix  A  (A. 62),  the  complete  expression 
for  the  MRTD  using  the  Ratches  Model  is: 
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MRTD 


—  7  NETD 
8 


■,Jl,HT(fy)Ho<fy)dfy 


vxAyi 


OS 


7 — —  HELECT  ( f x)Hd( f x.  f f X)HT( f y^C fy)d f xd f y 


(II. 8) 


To  obtain  this  equation  Ratches  et  al.  performed  the  following 
basic  steps: 

1.  The  actual  target  consisting  of  four  vertical  bars  is 
assumed  to  be  periodic  in  the  horizontal  direction. 

2.  The  human  eye-brain  system  is  represented  by  a  matched 
filter  consisting  of  one  bar  undegraded  in  the  horizontal 
direction  and  degraded  in  the  vertical  direction,  i.e., 

Hmf(fx.fy)=HT1bar(fx>fy)HD(fy>=HuCfx)HT(fy)Ho(fy)  (II.  9) 

3.  The  signal  is  obtained  by  taking  the  difference  between 
the  output  of  that  matched  filter  when  centered  over  the 
displayed  bars  and  when  centered  over  the  background. 
The  signal  is  referred  to  temperature  units. 

4 .  The  noise  in  the  display  is  assumed  to  be  white  and  it  is 
referred  to  temperature  units  by  means  of  the  NETD  and 
the  noise  equivalent  electrical  bandwidth  Afn. 

5.  The  signal-to-noise  ratio  present  at  the  display  is 
improved  by  two  factors  related  to  the  operator:  the  eye- 
brain  matched  filter  and  the  integration  performed  by  the 
eye  over  one  eye  integration  time. 

6.  Resolution  of  the  bars  is  considered  to  be  achieved  when 
the  signal-to  noise  ratio  has  exceeded  a  threshold  value 
Y.  In  the  case  of  Ratches'  paper,  the  experimental  value 
recommended  for  7  is  2.25  [Ref.  5],  [Ref.  8].  However, 
Lloyd  recommends  a  different  value  based  on  probability 
of  detection  of  single  bars  of  90  percent:  Y=4.5 
[Ref.  10,  p.  188].  Since  this  difference  does  not  affect 
the  shape  of  the  general  curve,  and  in  order  to  be 
consistent  with  the  plot  presented  by  the  latter,  the 
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simulation  of  the  Ratches  model  is  performed  with  the 
second  threshold  value  of  signal-to-noise  ratio. 

The  computer  simulation  of  (II. 8)  for  the  sample  system 
provides  the  plots  shown  in  Figure  II. 3  (linear  and  logarith¬ 
mic  scales) .  MRTDs  are  commonly  expressed  in  Celsius  degrees 
[Ref.  5]. 

C.  THE  LLOYD  APPROXIMATE  MODEL  FOR  THE  SAMPLE  SYSTEM 

The  authors  of  the  Night  Vision  Laboratory  model  [Ref.  5], 
present  some  approximations  to  (II. 8)  to  allow  hand  calcula¬ 
tions.  These  approximations  -acceptable  only  for  low  spatial 
frequencies-  are: 

LJlB5<VH»<VI^<Vdfy'1  '  (II. 10. a) 


LjlH"VI&VdV1<  and 


(Il.io.b) 


s 1  - 

2WJo  S7! - '  (II .  10  .  C) 

J  l(frx) 


since  from  definitions  of  Hy(fx)  and  HT(f  }  given  in  Appendix  C 


f_H?<Vdt»=X'  and 


(II. 11. a) 


P*S(fx>dfx=-^j 


(Il.ll.b) 


In  addition,  there  is  an  assumption  of  white  noise  at  the 
input  of  the  detector  preamplifier, 
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Figure  II. 3.  Predicted  MRTD  using  Ratches  Model 


(XX. 12) 


c  t 

'<v  _n 

e  I  ■L 

b  i^r*> 

Finally  the  response  of  the  system  in  the  low  frequencies 
satisfies: 

H0-Hd«HELECT«l  .  (11.13) 

By  using  these  approximations,  (II. 8)  reduces  to: 

MRn>=^THETDf  Ay.v,  2  fa  ll  (II.14) 

8  Hocw  [F^Tn  7  «j 

which  leads,  for  the  sample  system,  to  a  simple  approximate 
MRTD  expression: 

fox 

MRTD-0  •  49  _  ox  .  (11.15) 

Lloyd  [Ref.  10,  p.  190]  obtains  this  simplified  form  by 
working  strictly  in  the  temporal  frequency  domain.  This 
result,  (11.15),  convenient  for  some  applications,  is  less 
reliable  for  the  recognition  and  identification  task  at  higher 
spatial  frequencies  than  the  more  general  form  (II. 8).  The 
computer  simulation  of  (11.15)  is  presented  in  Figure  II. 4  and 
compared  with  Hatches'  predictions  (II. 8).  A  good  agreement 
is  observed  in  the  low  spatial  frequencies  while  the  effects 
of  the  approximations  become  noticeable  in  the  high  spatial 
frequencies. 

Note  that  according  to  (II. 15),  the  Ratches  model  [Ref.  5] 
predicts  that  the  MRTD  would  go  to  zero  as  fox  -♦  0.  Comparison 
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with  actual  measured  MRTD  data,  e.g.,  as  shown  in  Figure  1.9 
and  Figure  I. 10,  indicates  that  the  low  target  frequency  limit 
in  the  analytic  expression  (11.15)  is  not  representative  of 
actual  TIS  behavior.  This  is  one  of  the  motivating  factors 
for  the  development  of  the  visibility  model  proposed  in  the 
next  chapter. 


D.  THE  VORTMAN  AMD  BAR- LEV  MODEL  FOR  THE  SAMPLE  SYSTEM 

Vortman  and  Bar-Lev  [Ref.  8],  introduced  another  approach 
for  the  MRTD  model  based  on  Ratches.  The  main  difference  with 
respect  to  the  original  model  is  the  inclusion  of  an  adaptive 
matched  filter  to  represent  the  eye-brain  system.  This  filter 
is  matched  to  the  fully  degraded  target  signal  rather  than  the 
traditional  predetermined  matched  filter  presented  by  Ratches. 
Also,  in  the  Vortman  and  Bar-Lev  model,  the  target  is  consid¬ 
ered  in  its  actual  spatial  description  instead  of  assuming 
horizontal  periodicity. 

In  general,  any  model  based  upon  Ratches  that  uses  the 
actual  full  target  spectrum  rather  than  its  horizontally 
periodic  approximation  (A.  41)  has  an  MRTD  expression  where: 


8  4  wj>,,x)IWx)[e<'2''»“,-l]df, 


(II. 16. a) 


and,  therefore: 
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Figure  II. 4.  Predicted  MRTD.  Lloyd's  approximation  vs.  Ratches 
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YNETP 

at 


VxAy, 

TSI 


ELECT<fx>Hd<fx,fy)H«rf<fx,fy>dfxdfy 


,1 

2 


(II.16.b) 

Equation  (II.16.b)  is  a  generic  expression  for  the  MRTD 
compatible  with  the  Ratches  [Ref.  5],  Lloyd  [Ref.  10],  and 
Vortman  and  Bar-Lev  [Ref.  4]  models.  What  distinguishes  these 
models,  within  the  context  of  (II.16.b),  are  the  matched 
filter  (H^)  and  the  target  spectrum  (HT) . 

The  adaptive  matched  filter  employed  in  this  model  is  a 
single  bar  of  the  target  as  in  the  Ratches  model.  However,  it 
is  different  from  the  original,  only  degraded  in  the  vertical 
direction  (A. 60),  since  both  horizontal  and  vertical  degrada¬ 
tions  are  now  included.  Therefore,  in  the  spatial  frequency 
domain: 

H»f(fXffy)=HU(fx)HT(fy)Hopt(fXffy)Hdet<fx,fy)HELECT(fx)Hd(fx,fy)  (II.  17. a) 

Although  Vortman  and  Bar-Lev  include  a  transfer  function 
for  the  eye  for  (11.16),  this  is  not  considered  here  in  order 
to  maintain  consistency  with  the  other  models  analyzed. 

These  authors  also  discuss  the  consequences  of  employing 
white  noise  assumptions  in  the  model.  This  implies  working 
with  a  matched  filter  for  white  noise,  instead  of  a  theoret¬ 
ically  more  realistic  matched  filter  for  colored  noise.  The 
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Vortman  and  Bar-Lev  analysis  for  the  colored  noise  MRTD  led  to 
physically  untenable  conclusions,  so  the  model  was  discarded. 

Finally,  after  using  the  complete  spectral  description  of 
the  target  and  the  adaptive  matched  filter,  (II. 8)  is  trans¬ 
formed  into  Vortman  and  Bar-Lev  model  for  white  noise: 


MRTD 


TNETD 

at 


2‘AfnFrteJ-»J-«HELECT<fx>Hd{fK<fy>Hwffx>HT<fy>HoPt(f*-fy)Hdet(fK*fy)dfxdf' 


.  1 
2 


<,>"»<  VH'<V^<<x.V [  ]  df.df. 


(II.17.b) 

which  is  plotted  in  Figure  II. 5  versus  the  Ratches  model 
(II.  8).  It  can  be  observed  that  here  is  a  substantial 
disagreement  between  both  models,  specially  as  spatial 
frequency  increases. 

The  reader  can  confirm  that  (II.17.b)  is  in  agreement  with 
equation  16  from  [Ref.  8]  after  noting  that  td=Ax/vx. 

The  three  models  considered  in  this  chapter  represent 
different  concepts  of  visual  perception.  In  all  three  cases, 
discrepancies  between  predicted  and  measured  data  persist. 
The  need  for  improved  modeling  is  the  motivation  for  the 
proposed  models  introduced  in  the  next  chapter. 
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III.  PROPOSED  MRTD  MODELS 


A.  MODIFIED  VORTMAN  AMD  BAR-LEV  MODEL  FOR  THE  SAMPLE  SYSTEM 

As  presented  in  Appendix  A,  the  Ratches  model  is  based  on 
the  assumption  of  white  noise  at  the  input  of  the  TIS.  The 
transfer  function  of  the  system  obviously  colors  the  displayed 
noise  as  stated  by  Vortman  and  Bar-Lev,  and  therefore  the  use 
of  a  matched  filter  for  white  noise  (impulse  response  equal  to 
the  spatially  reversed  displayed  signal)  should  not  be 
considered.  However,  practical  reasons  require  the  assumption 
of  a  matched  filter  for  the  white  noise  case. [Ref .  8] 

If  the  input  noise  is  assumed  to  be  white  as  well  as  the 
displayed  noise,  this  implies  the  system's  MTF  does  not  affect 
the  noise  spectral  distribution.  It  follows  that  the  same 
criterion  must  be  employed  in  the  definition  of  the  matched 
filter,  and  no  degradation  should  be  considered  for  the 
original  bar  in  the  direction  where  the  noise  is  assumed  to 
occur  (horizontal) .  Therefore,  presenting  a  horizontally 
degraded  adaptive  matched  filter  under  white  noise  assumptions 
is  a  contradiction  that  leads  to  inaccurate  modeling.  This 
was  evident  in  the  computer  simulations  shown  in  Figure  1.7 
and  Figure  1.8. 

Two  different  modifications  can  be  attempted  on  the 
Vortman  and  Bar-Lev  model: 
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1.  Matched  Filter  undegraded  only  in  the  horizontal  direc 
tion  (where  the  considered  noise  occurs) : 

Hmf(f*,fy)=H«<f*)HT(fy)Ho(fy)  (Ratches)  (Ill.l.a) 


which  leads  to  a  modified  Vortman  and  Bar-Lev  MRTD  (Mod. 
No.l)  after  substitution  into  (II.16.b): 


MRTD  = 


7NETD 


2AfnFPteJ-® J1Helect<  fx)Hd(  fx,  f  y)H«(  f  X)HT(  f  y)H0(f  y)d  f  xdf  y 


JlJlHT<fx.VHu<fx>HTCVH0(fx.VH0(fy)[e  xU)-l]dfxdf, 


(Ill.l.b) 

where,  as  explained  in  the  early  part  of  Appendix  Af  the 
single-bar  normalized  spectrum  Hy{f  )  can  supplant  the 
four-bar  normalized  spectrum  HT(f  }  without  loss  in  gen¬ 
erality.  x 

2.  Matched  Filter  undegraded  in  both  horizontal  and  vertical 
directions : 


Hnrf<fx.V=HM<fX>HT(fy>  (III.  2.  a) 

which  leads  to  a  second  modified  Vortman  and  Bar-Lev  MRTD 
(Mod.  No.  2): 


MRTD  « 


YNETD 


2^f^^JT.JT»HELECT(fx)Hd(fx,fy>^<fx)^(fy)dfxd^ 


(III.2.b) 

Figure  III.l  shows  the  simulation  of  the  first  modified 
model.  The  predicted  MRTD  is  coincident  with  the  one  obtained 
by  Ratches. 
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MRTD 


Modif.  V-  B-L  No.1("),Ratches(-)  (Lin  Scale)  Hod  if.  V-B-L  No.1(~),Ratches(-)  (Log  Scale) 


Spatial  Frequency  (cy/arad)  Spatial  Frequency  (cy/arad) 


Figure  III.l.  Predicted  MRTD  using  modification  No.  1  on  V-B-L 
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Figure  III. 2  shows  the  second  modification,  where  the 
results  observed  are  essentially  the  same  as  for  the  first. 
This  demonstrates  that  the  issue  of  degradation  along  the 
vertical  direction  is  irrelevant  for  a  typical  TIS  like  the 
one  employed  in  this  work  and  that  the  Ratches  simplifications 
are  correct  in  that  sense. 

B.  A  VISIBILITY  MODEL  FOR  MRTD  MEASUREMENTS 

Today  fairly  successful  MRTD  predictors  appearing  in  the 
literature  [Ref.  5],  [Ref.  8],  [Ref.  9],  [Ref.  10,  p.  190]  are 
all  based  on  a  signal-to-noise  (S/N)  analysis.  The  predicted 
MRTD  is  subsequently  determined  by  a  threshold  condition  on 
the  S/N.  In  most  cases  these  models  generate  an  MRTD  which  is 
directly  proportional  to  the  NETD.  This  parameter  is  defined 
as  the  temperature  difference  between  a  large  target  and 
background  which  leads  to  a  S/N=l  under  conditions  that  a 
special  reference  filter  is  employed.  The  reference  filter 
has  the  effect  of  desensitizing  the  NETD  to  spatial  frequency 
effects  inherent  to  the  system.  Excluding  the  NETD,  the  rest 
of  the  MRTD  will  depend  on  inherent  spatial  frequency  effects. 
This  motivates  a  proposed  generic  form  for  the  MRTD: 

MRTDtfox)  =  MRTDtWO)  function<«yst«m  NTFs) 

which  does  not  appear  to  be  100%  consistent  (see  11.15)  with 
previous  MRTD  predictors  discussed  in  the  literature.  None¬ 
theless  this  form  has  significant  physical  appeal. 
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MRTD 


Mod i f .  V-B-L  No.2(--),Ratches(0  (Lin  Scale)  Mod i f .  V-B-L  No.2(--),Ratches(-)  (Log  Scale) 


Spatial  Frequency  (cy/arad) 


Spatial  Frequency  (cy/arad) 


Figure  III. 2.  Predicted  MRTD  using  modification  No.  2  on  V-B-L 
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In  the  low  frequency  limit  the  MRTD  tends  to  be  more 
optimistic  than  the  NETD  due  to  eye-brain  matched  filtering 
[Ref.  5],  [Ref.  8].  In  these  established  models  the  complex 
[Ref.  14]  and  controversial  [Ref.  25]  task  of  treating  eye- 
brain  filtering  has  to  be  addressed. 

In  the  visibility  model  to  be  discussed,  the  complexity  of 
eye-brain  modeling  is  sidestepped  by  not  dealing  with  the  NETD 
directly.  The  low-frequency  MRTD,  which  serves  partially  as 
an  NETD  substitute,  also  incorporates  effects  of  the  MRTD 
decision  process.  For  the  standard  subjective  MRTD  this 
involves  eye-brain  filtering.  The  analysis  for  the  visibility 
model  concentrates  on  contrast  reduction  due  to  spatial 
frequency  limiting  factors.  Direct  consideration  of  noise  is 
obviated . 

This  proposed  model  leads  to  a  simpler  development  for  an 
MRTD  predictor  and  appears  to  provide  a  more  reliable  tool  for 
forecasting  the  performance  of  a  TIS. 

As  pointed  out  by  Lloyd  [Ref.  10,  p.  210],  the  modulation 
transfer  function  "correlates  with  recognition  and  identifica¬ 
tion  performance  in  noise-free  imagery".  The  degradation 
produced  by  the  system  to  a  rectangular  wave  in  the  horizontal 
direction  is  considered  and  the  amplitude  of  the  resulting 
wave  is  compared  to  a  critical  value  to  determine  the  recogni¬ 
tion  of  the  target.  Hypothetically,  this  threshold  for 
recognition  could  eventually  be  expressed  as  a  function  of  the 
NETD  and  other  system  parameters. 
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Figure  III. 3  summarizes  the  salient  features  of  the  model. 
The  overall  MTF  is  obtained  from  subsystem  effects.  The 
thermal  rectangular  wave  source  is  characterized  by  a  period 
l/fox  and  temperature  difference  AT.  The  response’s  amplitude 
ATS  is  guaranteed  to  be  less  than  AT.  The  numerical  curves 
discussed  in  the  next  paragraph  were  obtained  by  applying  the 
computation  algorithm  shown  in  Appendix  B  on  the  sample  system 
described  in  Appendix  D. 


l/fox 

t=>  C=j 


_^JUUl  at 

0 


MTF 


Figure  III. 3.  Modified  Visibility  Model  Concept 

The  upper  half  of  Figure  III. 4  shows  the  approximation  to 
a  rectangular  wave  of  spatial  frequency  equal  to,  e.g.,  0.7 
[cy/mrad]  employing  the  first  five  harmonics  of  its  Fourier 
series  expansion  (undegraded  target) .  The  lower  half  of  III. 4 
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shows  the  degraded  version  of  that  signal  after  passing 
through  the  TIS  (degraded  target) .  Besides  a  loss  of  the 
original  shape ,  its  amplitude  has  decreased. 

Upper  and  lower  portions  of  Figure  III. 5  show  the  indi¬ 
vidual  5  harmonics  in  the  undegraded  and  degraded  targets 
respectively  (the  latter  shows  the  four  highest  harmonics 
essentially  canceled) . 

Since  the  degraded  target  is  the  available  signal  for  the 
observer,  it  provides  distorted  information  on  what  is  the 
actual  difference  in  temperature  in  the  target.  If  this 
measured  temperature  difference  is  called  ATS,  then  a  param¬ 
eter  a  that  relates  the  measured  temperature  difference  to  the 
actual  one  (AT)  can  be  defined: 


For  the  sample  system,  this  parameter  a  is  plotted  versus 
spatial  frequencies  in  Figure  III. 6. 

It  is  reasonable  to  assume  that  there  is  a  minimum  or 
critical  measured  apparent  temperature  difference  (AT#C)  that 
an  observer  can  resolve.  By  definition  AT#e  corresponds  to  the 
condition  AT=MRTD.  Thus: 


MRTD= _ —  .  (III. 4) 

a 

Figure  III. 7  shows  the  resulting  predictions  compared  with 
the  Ratches  model  for  the  sample  system,  with  AT#c=0.23*C. 
This  critical  temperature  difference  produced  fairly  close 
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Figure  III. 4.  MTF  degradation  of  a  rectangular  wave 
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agreement  with  the  Ratches  model.  In  order  to  relate  the 
visibility  model  to  previous  descriptions,  ATSC  can  be  ex¬ 
pressed  as: 

ATse=M,-P  (III.  5) 

0 

where  0  should  be  an  empirical  factor  dependent  on  the  mea¬ 
surement  process  for  recognition.  Then,  ATSC  clearly  becomes 
a  new  NETD-like  parameter  where  the  signal-to-noise  ratio 
improvement  performed  by  the  eye-brain  system  and  the  display 
effects  are  included.  For  the  simulation  of  this  model  on  the 
sample  system,  the  0  calculated  from  the  NETD  (see  end  of 
Chapter  II,  Section  A)  was  0=1.58. 

The  curves  shown  in  Figure  III. 7  demonstrate  several 
noteworthy  trends.  At  low  frequencies  the  visibility  model  is 
not  as  optimistic  as  the  Ratches  model.  At  higher  frequencies 
this  relation  is  reversed  and  the  visibility  model  tends  to  be 
more  optimistic  than  the  Ratches  model.  Figure  1.9  and  Figure 
I. 10  show  that  measured  data  behaves  similarly  when  compared 
with  the  Ratches  model.  Therefore  the  visibility  model  has 
excellent  potential  for  serving  as  a  more  reliable  performance 
forecaster  than  the  Ratches  model. 

C.  INCLUSION  OF  SAMPLING  ARTIFACTS.  EFFECTS  ON  THE  RATCHES 

MODEL 

For  systems  considered  here,  as  represented  in  Figure  A.l, 
the  main  sampling  effect  takes  place  along  the  vertical 


41 


MRTD 


Visibility  Mod.(--)»Ratch«s(-)  (Un  Scale)  Visibility  Kod.(— )JRatches(-)  (Log  Scale) 


Spatial  Frequency  (cy/arad)  Spatial  Frequency  (cy/arad) 


Figure  III. 7.  Visibility  MRTD  model 
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direction.  This  is  due  to  the  horizontal  scanning  being 
performed  on  lines  separated  by  an  angular  distance  mrad. 
This  consequently  leads  to  a  periodicity  of  the  spectrum  of 
the  observed  target  in  the  vertical  direction.  The  target 
spectrum  repeats  itself  every  1/Ay,.  [cy/mrad].  In  Appendix  C, 
the  full  2D  spectrum  of  the  four  bars  is  derived.  It  is  noted 
that,  consistent  with  specifications  in  Appendix  D  for  the 
sample  system,  the  center-to-center  spacing  of  the  detectors 
is  taken  to  be  equal  to  the  vertical  detector  width.  In  a 
well-designed  system,  the  main  energy  carried  by  the  aliases 
crated  by  sampling  should  lay  in  the  region  beyond  the  cutoff 
frequency  of  the  overall  MTF.  In  order  to  evaluate  the 
importance  of  this  unwanted  effect,  a  periodic  spectrum  is 
considered  here  for  MRTD  modelling.  Figure  III. 8  shows  the 
periodic  (dashed)  and  aperiodic  (solid)  vertical  spectrum  of 
a  standard  target  of  0.7  [cy/mrad)  for  the  sample  system 
(Appendix  D) .  Figure  III. 9  shows  how  the  system's  MTF 
minimizes  this  effect. 

As  an  example,  Figure  III. 10  shows  the  plot  of  the  MRTD 
according  to  Ratches,  using  a  the  new  sampled  HT(fy).  In  this 
case,  as  expected,  there  is  no  essential  difference  with 
respect  to  the  original  model.  This  shows  that  the  sample 
system  is  correctly  designed  since  the  modulation  transfer 
function  is  such  that  the  periodic  spectra  are  essentially 
filtered  out  and,  therefore,  sampling  effects  can  be  neglect¬ 
ed. 
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Figure  III. 9.  MTF  effects  on  vertical  spectrum 


MRTD 


Batches  t/sanpl .(--), Ratches(-)  CUn  Scale)  Ratches  */sanpl Ratches(-)  (Log  Scale) 


Spatial  Frequency  (cy/»rad)  Spatial  Frequency  (cy/«rad) 

Figure  III. 10.  Sampling  artifacts  on  Ratches'  MRTD  prediction 
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IV.  COMPARATIVE  ANALYSIS  AND  CONCLUSIONS 

Some  of  the  conclusions  presented  here  are  based  on 
computer  simulations  performed  on  the  sample  system  specified 
in  Appendix  D. 

As  expected,  the  Lloyd  model  demonstrates  the  best 
agreement  with  the  Ratches  model  at  low  spatial  target 
frequencies.  However,  as  observed,  there  is  a  physical 
inconsistency  with  the  Ratches  model  at  low  frequencies.  This 
led  to  the  generally  recognized  conclusion  that  a  model  that 
provides  a  higher  estimation  for  the  MRTD  at  low  frequencies 
is  needed. 

The  Adaptive  Matched  Filter  (AMF)  model  by  Vortman  and 
Bar- Lev  employs  a  mathematically  more  technical  application  of 
the  matched  filter  treatment  for  the  eye-brain  than  does  the 
Ratches  model.  Comparison  between  simulations  and  actual 
measured  data  suggest  that  Vortman  and  Bar-Lev  treatment 
produces  an  unwanted  and  magnified  disagreement  between  model 
and  reality  in  the  high  frequencies  without  obtaining  any 
improvement  in  the  lower  ones.  The  conclusion  from  this 
comparison  is  that  the  eye-brain  is  better  modeled  as  a 
matched  filter  for  the  undegraded  target  signal,  as  originally 
proposed  30  years  ago  by  Ratches. 

It  was  found  in  Chapter  II  and  Chapter  III,  by  comparison 
of  computer  simulations  for  the  sample  system,  that  use  of  the 
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exact  four-bar  target  spectrum  is  not  really  necessary.  Not 
only  can  the  spectrum  be  reasonably  approximated  by  an 
infinite  square  wave  but  Hatches'  use  of  only  the  first 
harmonic  turns  out  to  be  quite  good.  This  was  demonstrated  in 
Chapter  III,  Section  A,  by  redefining  Vortman  and  Bar-Lev 
matched  filter  (both  modifications  No.  1  and  No.  2)  .  The 
results  from  the  new  models  essentially  coincide  with  the 
Ratches  model.  This  proves  that  the  periodic  assumption  for 
the  target  and  the  use  of  the  first  harmonic  proposed  by 
Ratches  (A. 41)  is  a  correct  engineering  approach.  The  results 
obtained  from  modification  No.  2  to  Vortman  and  Bar-Lev  con¬ 
firm  the  general  irrelevancy  of  including  degradation  effects 
along  the  vertical  axis  within  the  matched  filter  as  in  the 
Ratches  model. 

The  proposed  Modified  Visibility  model  brings  up  a  very 
simplistic  concept  of  the  perception  process.  Although  the 
feasibility  of  a  universally  valid  parameter  fi  requires 
additional  research,  the  results  seem  to  be  in  better  agree¬ 
ment  with  physical  measurement  than  the  Ratches  model.  Since 
this  has  been  only  tested  on  the  sample  system,  further 
experimental  simulations  on  different  systems  must  be  pursued. 
It  should  be  noted  that  this  model  does  not  suffer  from  the 
unphysical  condition  that  the  low  frequency  limit  for  the  MRTD 
equals  zero.  It  was  demonstrated  in  Chapter  III,  Section  C 
that  the  Ratches  model  exhibits  this  defect.  The  possibility 
that  this  problem  is  endemic  to  all  models  based  on  a  S/N 
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calculation  needs  to  be  studied.  Furthermore  the  general 
development  of  the  visibility  model  makes  no  assumptions  on 
the  nature  of  the  decision  process.  For  this  reason  the  modi¬ 
fied  visibility  model  has  significant  unexplored  potential  as 
a  mathematical  framework  for  objective  MRTD  modeling. 

Finally,  it  was  demonstrated  in  Chapter  III,  Section  C 
that  the  inclusion  of  sampling  artifacts,  tested  on  the  sample 
system,  hardly  affected  the  MRTD  calculation.  In  order  to 
evaluate  the  impact  sampling  effects  will  have  on  the  MRTD 
calculation  for  an  arbitrary  system,  the  approach  presented 
in  Chapter  III,  Section  C  should  be  applied. 

For  readers  interested  in  a  more  detailed  description  of 
the  computer  simulation  algorithm  a  flowchart  is  provided  in 
Appendix  B. 
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APPENDIX  A  -  INTERPRETATION  OF  THE  RATCHES  MODEL 


The  inclusion  of  this  appendix  serves  two  important 
functions.  First,  it  provides  the  necessary  background  to 
make  this  thesis  fairly  self-contained.  Second,  the  evolution 
of  concepts  surveyed  here  includes  additional  physical 
interpretations  not  provided  in  the  original  manuscript 
[Ref.  5].  From  a  pedagogical  point  of  view  these  additional 
enhancements  are  essential  to  the  fundamental  understanding  of 
the  original  work. 

The  model  presented  by  the  NIGHT  VISION  LABORATORY  in  1975 
[Ref.  5],  better  known  as  the  Ratches  model,  predicts  system 
Noise  Equivalent  Temperature  Difference  (NETD) ,  Modulation 
Transfer  Function  (MTF) ,  Minimum  Detectable  Temperature 
Difference  (MDTD) ,  and  Minimum  Resolvable  Temperature  Dif¬ 
ference  (MRTD)  for  any  scanning  thermal  imaging  system.  An 
important  characteristic  of  this  model  is  the  fact  that  it 
represents  the  image  processing  that  takes  place  in  the  eye- 
brain  system  by  a  matched  filter  [Ref.  5],  [Ref.  8].  Although 
this  assumption  may  not  exactly  correspond  to  what  really 
occurs  in  the  complex  process  of  visual  recognition,  [Ref.  12], 
[Ref.  13],  [Ref.  14],  [Ref.  15],  it  has  shown  a  remarkable 
success  and  has  been  widely  accepted  by  the  infrared  communi¬ 
ty.  The  present  appendix  follows  the  steps  presented  by 
Ratches  and  Lawson  in  the  report  "The  Fundamentals  of  Thermal 
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Imaging  Systems"  produced  by  the  NAVAL  RESEARCH  LABORATORY  in 
1979  [Ref.  30]. 

A.  PRELIMINARIES 

From  linear  systems  theory,  it  is  known  that  the  output 
signal  from  a  system  is  egual  to  the  input  signal  convolved 
with  the  impulse  response  of  the  system,  i.e., 

i0{t)  =  *  k<t>  =  ^Cr>  h(t-r)CiT*  (A.l) 

where  the  subscript  i  stands  for  input  and  the  subscript  o  for 
output . 

In  the  freguency  domain,  the  Convolution  Theorem  can  be 
applied  to  (A.l)  to  obtain  the  expression: 

Io(f)  =  Ii(f)  H(f)»  (A. 2) 

where  I0(f)  is  the  Fourier  transform  of  i0(t)  and  H(f)  is  the 
transfer  function  of  the  linear  system. 

If  cascaded  systems  are  considered,  the  Fourier  transform 
of  the  output  signal  is: 

I0(f)  *  Ij<f)  Huo  Hgcf)***  (A. 3) 

where  H1(f)  and  H2(f)  are  the  transfer  Functions  of  the  first 
and  second  systems  respectively. 

To  characterize  the  effects  of  noise  generated  by  a  wide- 
sense  stationary  process  such  as  that  found  in  most  thermal 
imaging  systems,  the  autocorrelation  function  of  the  noise 


51 


needs  to  be  employed.  This  operation,  for  finite  power 
signals,  is  defined  as: 

T 

R,(„=u%.fV;„fOT*  <».«.»> 

J_2 

which  can  be  represented  for  the  real  noise  signal  n(t) 

R(r)  =  <  n(t>  >  •  (A.4.b) 

where  brackets  indicate  average.  A  random  process  is  said  to 
be  wide-sense  stationary  if  the  statistics  corresponding  to 
average  and  autocorrelation  do  not  depend  on  the  location  of 
the  time  origin  [Ref.  31].  Expressions  (A. 4)  have  already 
incorporated  this  assumption. 

The  power  spectral  density  (PSD)  of  the  noise  can  then  be 
found  by  taking  the  Fourier  transform  on  the  autocorrelation 
function  (Wiener-Khintchine  Theorem  [Ref.  31,  p.  458]): 

S(f)  =  9  {R(t)).  (A. 5) 

After  passing  through  a  linear  system,  the  output  PSD  is 
[Ref.  32]: 


=  Si(f) 


(A.  6) 


The  time  variance  of  a  measured  signal  due  to  the  noise 
generated  by  a  random  process  of  mean  equal  to  zero  is  related 
to  the  noise  spectral  density  as  follows: 
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(A. 7. a) 


or,  in  the  frequency  domain, 


(A.7.b) 


To  avoid  the  use  of  negative  frequencies,  a  one-sided  noise 
PSD,  S'(f),  is  defined  such  that: 


a2  =  jV<f>  df.  (A. 8) 

A  matched  filter  [Ref.  30]  "is  a  filter  whose  response 
function  is  a  delayed  (shifted) ,  time-reversed  (spatially 
reversed)  version  of  the  signal".  Thus,  if  the  signal  is 
defined  as  i(t),  the  impulse  response  of  a  matched  filter 
results: 

hm(t)  *  K  i(r  -  t>  (A.  9) 

where  K  is  a  constant.  For  convenience  this  arbitrary  con¬ 
stant  will  be  taken  to  have  the  ideal  value  1  and  units  of 
(time)*1.  This  kind  of  filter  maximizes  the  signal-to-noise 
ratio  at  time  t=r,  considering,  a)  the  signal  output  is  the 
peak  of  the  autocorrelation  of  the  input  signal  and  b)  the 
input  noise  that  assumed  to  be  additive  and  white.  This 
assumption  is  approximately  true  for  thermal,  shot,  genera¬ 
tion-recombination  and  radiation  (or  photon  flux)  types 
of  noise  [Ref.  33]. 
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Any  1/f  noise  component  in  the  low  frequencies  and  the 
high-frequency  roll-off  are  ignored  in  the  model.  Assuming 
r=0,  the  transfer  function  of  such  a  matched  filter  is: 

«W«-KICV  <A*10> 

The  final  MRTD  expression  depends  only  on  the  magnitude  of  the 
matched  filter.  Since  linear  transform  theory  predicts  that 
the  magnitude  of  a  transfer  function  does  not  depend  on  a 
shift  in  the  corresponding  impulse  response,  the  assumption  is 
justified. 

B.  DERIVATION  OF  THE  NOISE  EQUIVALENT  TEMPERATURE  DIFFERENCE 
(NETD) 

The  definition  of  NETD  employed  in  this  model  is  [Ref.  30] 
"that  input  temperature  difference  for  a  ’large’  target  (a 
large  target  being  one  whose  size  is  large  relative  to  the 
system  response  function)  which  is  required  to  generate  a 
signal  (voltage  amplitude)  just  prior  to  the  display  (or  after 
the  detector  preamplifier)  which  is  just  equal  to  the  rms 
noise  (voltage)  at  that  point,  assuming  that  the  filtering 
action  of  the  electronics  prior  to  the  measurement  point 
corresponds  to  that  of  a  'standard'  filter."  In  other  words, 
the  NETD  is  the  temperature  difference  between  a  large  uniform 
target  (which  implies  low  spatial  frequency)  and  a  uniform 
background  that  produces  a  signal-to-noise  ratio  equal  to  1  at 
that  point. 
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It  can  be  assumed  (Figure  A.l)  that  each  detector  plus  the 
associated  preamplifier  form  a  linear  system  with  a  impulse 
response  function: 

^OET-€lECT(A,t)  =  ra,t>  *  ^ELECT(t)  [ VOlts/watt ]  ,  (A.  11) 

where  r(A  t)  is  the  response  of  the  detector  [volts/watt]  and 
hELECT(t>  represents  the  response  of  the  preamplifier  and  the 
associated  circuitry.  The  signal  incident  onto  the  detector 
can  be  represented  as  a  function  of  X  and  t,  in  which  the 
wavelength  and  time  dependence  are  separable: 

pu,t>  =  i(t>  [watts/ Mm]  (A.  12) 

The  quantity  P(A  t)  can  be  considered  as  the  time-varying 
spectral  radiant  flux  falling  on  a  detector.  The  function  h<Plx) 
is  the  wavelength  dependent  part  [watts/ Mm] ,  and  i(t)  a 
normalized  (peak  value  equal  to  one)  time  function.  The 
response  at  the  detector  is  then  obtained  by  extending  the 
linear  systems  result  (A. 11)  and  integrating  over  wavelength. 
The  output  voltage  v#(t)  from  each  detector-amplifier  system  is 
obtained  as: 

Vs(t)  =  *  r(A,t)  *  ^ELECT(t)^^  [volts]  (A.  13) 

which  corresponds  to  a  detector-amplifier  system  like  the  one 
presented  in  Figure  A.l. 

To  simplify  the  expression,  it  can  be  assumed  that  the 
one-dimensional  transform  of  r<A  t)  is  separable  into  a  wave¬ 
length  and  a  frequency  dependent  part  R(Jk)  and  R(f),  such  that: 
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Detectors  Amplifiers 


Figure  A.l.  Detector-Amplifier  system 


where : 


and 


In  expressions  (A. 14),  fr  is  a  convenient  reference  frequency. 
A  specific  guideline  for  choosing  this  reference  frequency 
will  be  discussed  later  in  this  Appendix. 
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Then,  taking  the  one  dimensional  Fourier  transform  of 


(A. 13)  and  inserting  (A. 14. a)  leads  to: 


V8(f)  =  ^0(A)I(f)K(>)R(f)HELECT(f)dA.  [volts/Hz )  (A.  15) 

which  follows  from  the  linearity  property.  Simplifying  (A. 15) 
produces : 

V.(f,  =  i<f )^(f )^£LECT(f jJ'q  ^(A)^CA)^  [volts/Hz].  (A. 16) 

Therefore 


vs(t)  =  1  ’ <t>J^  ^(A)R(A)d>.  [volts],  (A.17) 

where : 

i'(t)  =  i(t)  *  r(t)  *  ^ELECT(t)  •  (A.  18) 

and  where  r(t)  is  the  inverse  transform  of  R(f). 

Inherent  to  the  definition  of  NETD  is  the  concept  of 
instantaneous  signal-to-noise  ratio  which  is  represented  as: 


S  _  vs<t> 
N  a 


(A. 19) 


where  a  is  the  rms  noise  voltage  (A. 7)  and  v#(t)  is  the  cor¬ 
responding  signal  after  the  preamplifier.  From  (A. 6)  and 


(A. 8) : 


1 

a=[J^  S  '  l(f)HELECT(f)^f  1  ^ 


(A. 20) 


and  therefore, 
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(A. 21) 


S  i'(t>^a>Ra>^ 

N  T* 

Equation  (A. 21)  provides  the  basis  for  the  NETD  calculation  by 
setting  the  signal-to-noise  ratio  equal  to  one.  However,  it 
needs  to  be  recast  as  a  function  of  more  standard  parameters, 
including  (obviously)  the  temperature  difference  between  the 
target  and  background  AT. 

The  authors  of  the  model  included  here  a  very  reasonable 
assumption  by  taking  i'(t)=l.  This  is  qualitatively  explained 
by  the  fact  that  the  normalized  time  function  i(t)  must  be 
equal  to  one  at  the  midpoint  of  the  signal  (scanned  target) 
where  i'(t)  is  measured.  As  represented  on  Figure  A. 2,  for  the 
horizontal  scanning  of  the  target,  the  duration  of  i{t)  is  much 
greater  than  the  response  functions  r(t)  and  hELECTCt)-  Therefore 
r(t)  and  hELECT(t)  are  assumed  to  have  approximately  delta- 
function  like  behavior,  which  leads  to  the  condition  i'(t)=l. 

In  order  to  recast  (A. 21),  the  concept  of  Spectral  Noise 
Equivalent  Power  (NEP{X))  must  be  introduced: 

NEP»>  -  ^  (».«) 

a 


where : 
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Figure  A. 2.  Sample  i(t),  r(t),  and  hELECT(t) 

vsa)  =  ^(a)r<a)  [volts/ nm] .  (A. 23) 

Note  that  since  the  NEP  is  defined  for  constant  optical  flux, 
the  time  dependence  in  this  context  has  been  dropped. 
Equation  (A. 22)  implies  a  linear  relation  between  the  optical 
power  and  the  resulting  voltage  signal,  as  shown  in  Figure 
A. 3.  The  slope,  R(A),  of  this  linear  relation  depends  on  wave¬ 
length.  Definition  (A. 22)  is  in  agreement  with  the  more 
standard  NEP  definition  in  which  a  narrow  band  filter,  AX,  is 
employed  [Ref.  34]. 
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of  the  circuitry  considered  in  the  noise  measurement  and  Ad  is 
the  orea  of  the  detector  [cm2] .  If  HELECT{f)  represents  the 
overall  transfer  function  of  the  preamplifier-measurement 
device  system,  the  noise  bandwidth  Afn  is  defined  by 


Afn  -  r|^HfLECI(„df 

Jot>  i(fr) 


(A. 25) 


such  that  by  direct  comparison  with  (A. 20): 

a=S'j(fr,Afn.  (A.  26) 

Therefore,  after  substituting  (A. 22),  (A. 23),  and  (A. 26)  into 
(A. 24),  it  is  found  that: 


A  2R 
Ad  K(x> 


(*) 


(A. 27) 


<S*iCfr>)  " 

This  leads  to  a  new  expression  for  the  signal-to-noise  ratio: 


S 

N 


^<a)d;x)<u 


J^A0(A,D(A)dX 


C*d  (7  sr^^LECKf ,df ]  2  [ A^AfJ 

Jot*  i(fr) 


(A. 28) 


Next,  the  model  assumes  for  a  simple  imaging  system: 


^(x)  = 


(A. 29 ) 


where : 

rj0(A)  =  Optical  efficiency  of  the  viewer, 

F  =  f/number=  focal  length  (d) /diameter  of  the 

aperture  stop  (D)  [Ref.  33,  p.  179], 

T  =  Temperature, 

L(a)  =  Spectral  Radiance  [watts/ (cm2  steradian  /xm)  ]  from 
the  source. 
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This  follows  after  noting  that  the  approximate  solid  angle 
subtended  by  the  lens  {with  respect  to  its  focal  point)  is: 


The  expression  (A. 29)  describes  the  focusing  process  performed 
by  the  lens  and  the  subsequent  collection  of  the  energy 
focused  onto  the  focal  plane  by  each  detector.  It  also  implies 
that  the  spectral  irradiance  onto  the  detector  H(i)  is: 

H<A)=“^  =  ^r?0<A)“^AT  twatts/(cm2 /im)  ] .  (A. 30) 


Finally,  after  substituting  (A. 29)  into  (A. 28): 

1 

s  vl ATf^°<>> D<*>dX 

_  -  . 

4F2(Afn)  2 


(A. 31) 


By  setting  the  signal-to-noise  ratio  equal  to  one,  AT 
becomes  the  NETD: 

i 

4  F2 ( Af  )  2 

NETD  =  - - - : *- -  [*C].  (A. 32) 

There  is  one  additional  condition.  By  convention  the  NETD 
parameter  assumes  the  use  of  an  external  measurement  filter. 
In  essence  this  restricts  the  measured  signal-to-noise  ratio 
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to  be  dependent  on  parameters  such  as  the  detector  angular 
subtense  (DAS) ,  but  not  the  bandlimiting  components  making  up 
the  TIS.  The  external  filter  should  be  designed  such  that  the 
overall  MTF  will  satisfy: 


H 


2 

ELECT(f)- 


1 

l+(^)2 


(A. 33. a) 


where  the  3-dB  power  bandwidth  meets  the  condition: 


where  rd 

nition 

that: 


f„=. 


2T, 


(A.33.b) 


is  the  detector  dwell  time.  It  follows  from  the  defi- 
(A.  25)  and  the  white  noise  condition,  S- '  >/Sf  •  tfr3=l , 


Afn=frl  (A.33.C) 

which  was  determined  through  use  of  a  convenient  table  of 
integrals  [Ref.  35].  For  consistent  comparison  of  NETD 
associated  with  thermal  imaging  systems  the  noise  bandwidth 
should  satisfy  (A.33.c). 

For  a  single  detector  system,  with  100%  scanning  efficien¬ 
cy,  the  dwell  time  is  given  by: 


r  =  1  Ax  Ay 

d  Fr  (HFOV) (VFOV) 


(A. 34 ) 


which  follows  from  the  geometry  shown  in  Figure  1.4.  Ex¬ 
tending  this  result  to  cover  non-unity  scanning  efficiency  and 
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overscan  ratio  as  well  as  systems  based  on  a  parallel  detec 
tion  scheme,  leads  to: 


1  _  (HFOV)  (VFOV)  Frr?ovsc 
r~2Td  2npAxAyr?sc 


(A. 35) 


where : 

HFOV  =  device  horizontal  field  of  view  [mrad] 

VFOV  =  device  vertical  field  of  view  [mrad] 

Fr  =  frame  rate  [Hz] 

T70vsc  =  overscan  ratio  for  the  device 
np  =  number  of  detectors  in  parallel 
Ax  =  horizontal  detector  size  [mrad] 

Ay  =  vertical  detector  size  [mrad] 

tj8C  =  scan  efficiency  (fraction  of  time  spent  in  actually 
scanning  the  field) ,  and 

rd  =  dwell  time  (time  the  detector  takes  to  scan  a 
picture  element  AxAy)  [sec] . 

Again,  it  is  important  to  remember  that  the  actual  system 
noise  bandwidth  is  the  one  stated  in  (A. 25),  and  that  the 
approximation  given  in  (A.33.C)  is  only  valid  if  the  measuring 
device,  which  includes  H2ELECT(f),  is  adjusted  so  that  the  true 
bandwidth  of  (A. 25)  equals  the  standardized  one  of  (A.33.C). 
A  demonstrative  plot  of  Afn  and  H2ELECT{f)  for  the  sample  system 
employed  in  this  thesis  is  presented  in  Figure  A.  4.  The 
physical  interpretation  is  that  the  area  under  the  curve  is 
equal  to  Afn. 
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Afn  t=C> 

frequency  (KHz) 

Figure  A.  4.  Noise  equivalent  bandwidth  (Afn)  and  H2ELECT(f) 

A  standard  approximation  which  is  usually  used  to  simplify 
(A. 32)  is  introduced  in  Chapter  II,  Section  A. 

C.  DERIVATION  OF  THE  MINIMUM  RESOLVABLE  TEMPERATURE  DIFFER¬ 
ENCE  (MRTD) 

The  standard  target  with  which  the  MRTD  has  historically 
been  measured  and  modeled  consists  of  four  bars  in  which  the 
height  (L)  is  equal  to  seven  times  their  width  (W) ,  equally 
spaced  by  a  distance  W.  The  derivation  for  the  MRTD  that 
follows  considers  this  standard  pattern  and  assumes  no 
sampling  effects  along  the  direction  of  scanning.  Probably 
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the  most  important  assumption  of  this  derivation  is  the  repre¬ 
sentation  of  the  perceptual  recognition  process  that  takes 
place  in  the  eye-brain  system  by  means  of  a  matched  filter. 

For  a  viewed  object  (the  image  obtained  from  the  system) 
characterized  by  the  dimensionless  spatial  function  i(xy),  such 
a  matched  filter  would  have  a  response  function: 

kfliffx.y)  =  ^  ^(-x,-y)  (A.  36) 

where  K  is  a  constant  that  can  for  convenience  be  set  equal  to 
one,  but  will  carry  the  standard  units  for  a  2D  space  domain 
matched  filter  [angle'2] . 

From  the  space  domain  version  of  (A. 6),  the  output  rms 
noise  from  that  matched  filter  would  be: 


N*f  "  Cf^S(fx.V^f(fx.fy)dfxdfy]'2  <*.37.a) 

where  is  the  Fourier  transform  of  the  matched  filter 

impulse  response  (A. 36).  It  follows  that: 


Nf 


(A.37.b) 


where  S,.  .  .  is  now  the  double-sided  PSD  of  the  observed  noise 

lTX#TyJ 

and  I (fx>fy)  is  the  2D  spatial  Fourier  transform  of  i(xy). 

Mainly  because  of  the  MTF  degradations  that  must  be 
considered,  the  authors  explain  the  computational  advantage  of 
assuming  the  matched  filter  as  one  established  for  a  " (poten¬ 
tially)  infinite  periodic  pattern". 
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In  the  direction  of  scanning  (x)  the  matched  filter  is 
defined  as  an  undegraded  rect  function  of  width  equal  to  the 
width  of  the  bars.  In  the  other  direction  (y)  it  is  just  the 
degraded  rect  function  of  width  equal  to  the  length  of  the 
bars.  This  can  be  explained  by  thinking  of  the  human  eye- 
brain  system  as  one  prepared  to  recognize  four  vertical  bars. 
The  horizontal  direction  being  critical,  the  observer  expects 
sharp  edges  along  x  (regardless  of  the  sharpness  in  y)  to 
accomplish  the  recognition  process.  Figure  A. 5  shows  this 
concept . 


Because  both  signal  and  noise  are  passed  through  the  same 
matched  filter  the  choice  to  use  a  single  bar  representation 
is  only  one  of  mathematical  convenience. 
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To  perform  the  MRTD  derivation,  Ratches  et  al.  started  by 
defining  the  signal  energy  per  unit  angular  area  of  the 
display  as: 


k  AT  i 


(x,y) 

Yi  vx 


,  Joules , 

l - =-  J » 

mrad* 


(A. 38) 


where : 

kAT  =  watts  emitted  by  a  display  element  for  a  large 

target  with  a  temperature  difference  AT  with  respect 
to  background, 

Ayi  =  angular  distance  between  scan  lines  {by/r)msc)  , 
vx  =  angular  scan  velocity  of  display  element  (A x/rd) , 
i<x  y)  =  sPat*-ai  distribution  of  viewed  object 

=  ^-Ttx.y)  *  ^(x.y)' 

t d  =  dwell  time, 

iT(x  y)  «  spatial  distribution  of  the  original  target,  and 
hp^  y)  =  system  response  function. 

Equation  (A. 38)  provides  a  useful  rule  for  conversion  between 
temperature  and  energy  units.  Proportionality  with  AT 
reflects  the  AC  coupling  inherent  in  most  FLIR  systems. 

Formally  expressed,  the  assumed  periodic  target  iT(xy) 
results : 

iT(x.y)  =  s9n(sin(2,rfoxx)  )  iuy)'  (A. 39) 

where  f0K  is  the  horizontal  spatial  frequency  of  the  target 
measured  in  cycles/mrad.  This  implies  that  the  function  iT(x  y) 
can  be  separated  into  two  independent  functions  along  each 
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Figure  A.  6  depicts  a 


dimension  such  that  iT(x,y)!=^T(x>^T<y)* 
sample  iT(x). 


Figure  A.  6.  Sample  iT(x) 

The  periodic  model  for  the  target  is  acceptable  if  the 
output  signal  is  seen  as  the  difference  between  the  energy  of 
the  target  passed  a)  through  the  matched  filter  centered  over 
the  target  and  b)  through  the  matched  filter  centered  over  the 
background . 

Applying  (A. 38),  and  the  concept  of  periodic  matched 
filter,  the  signal  used  in  the  MRTD  calculation  will  be: 
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S=Max|^-^i(x<y)*h^(Xjy)j-Min|^-— -i(*y)*h1frf(X>y) 

representing  the  dynamic  range  in  the  energy  signal. 

In  order  to  simplify  the  derivation,  the  model  assumes 
that  the  function  i(x>y)  can  be  approximated  horizontally  by 
using  the  first  two  terms  of  the  series  expansion  of  iT(x).  The 
resulting  sine  expansion  chosen  by  the  authors  corresponds  to 
an  odd  iT(x).  An  even  function  is  used  in  Appendix  c  which  is 
consistent  with  a  centered  target  in  the  object  plane.  This 
difference  is  not  relevant  for  the  final  result.  Thus: 


r  Joules i 
mrad2 


(A. 40) 


i(K.y)=t°*5  +  «0 <fox)  |  sin(2tff0Xx)  ]  i(y)  (A. 41) 

where  i(y)  is  the  degraded  vertical  function  that  corresponds 
to  the  length  of  the  target.  Therefore,  substituting  (A. 41) 
into  (A. 40)  leads  to: 


Joules j 
mrad2 


(A. 42 ) 


In  the  calculation  of  (A.  42)  the  DC  term  in  (A. 41)  can¬ 
celed  out.  Both  Min[]  and  Max[]  terms  in  (A. 40)  contribute 
equally.  As  discussed  in  Chapter  IV,  the  replacement  of  the 
full  periodic  spectrum  by  the  first  two  terms  turns  out  to  be 
an  extremely  good  approximation  for  the  sample  system  speci¬ 
fied  in  Appendix  D. 
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In  equation  (A. 42) ,  the  horizontal  integral  is  evaluated 
in  space  and  the  vertical  in  spatial  frequency.  The  first 
integral  uses  a  matched  filter  of  width  l/(2fox)  and  an 
amplitude  of  2fox  under  the  convention  of  a  normalized  area 
of  one  (H(f)i)=l  for  fx=0)  .  The  first  integral  equals  2/n  and 
the  second  one  can  be  simplified  since: 

Hmftfy)  =  HT(fy)  HflCfy)  (A.  43.  a) 

and 

■^(fy)  =  L  HT(fy)  HtXfy)  (A.43.b) 

where : 

L  =  length  of  a  bar  [mrad] , 

HT(fyj  =  normalized  transform  of  the  target  along  the  y 
direction  (see  (C.ll)),  and 

Hp^yj  =  normalized  transfer  function  of  the  overall  system 
along  the  y  direction. 

Finally  the  expression  for  signal  considered  for  MRTD  can 
be  recast  to: 


<*•«» 

With  respect  to  the  noise  analysis  that  will  be  used,  the 
model  starts  by  recognizing  the  need  to  produce  an  expression 
for  the  power  spectrum  of  the  displayed  noire.  The  noise  on 
the  display  is  given  by  the  following  equation: 
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n<x,y>-]CbHx)5<y-yj)*hd<y)  "  52  bi(x)bd<y-yj>  t  — - - 3  (A.  45) 

i  i  mrad* 

where  hd(y)  is  the  impulse  response  of  the  display  in  the  y 
direction.  The  function  bj(x)  is  related  to  a  voltage  noise 
function,  along  the  ith  video  line,  which  has  been  converted 
to  a  one-dimensional  radiant  function  by  the  display,  and  the 
spread  produced  on  each  line  of  the  monitor  (represented  by 
the  convolution)  .  Therefore,  n{x<y)  includes  the  noise  contri¬ 
bution  of  all  the  scanning  lines  on  every  single  point  of  the 
display.  The  autocorrelation  of  n(x  y)  is: 


<n<x,y)n(xT,yT)>=<52  52  bUx)bj(xT)hd(y-yi)hd(yT-yj)> 

i  1  Tnn1oe  (A. 46) 

=52  52  <b»<x)bj(xT)>bd(y-yi)bd(yT-yj)  t - ~2“  ]2 

•  •  J  mrad* 


i  j 


since  <[]>  represents  averaging  along  the  x  direction,  where 
the  random  process  takes  place. 

Ratches  et  al.  assume  a  noise  with  mean  equal  to  zero 
(<bi(x)>=0)  .  Then  the  result  of  <bj(x)  bj(XT>>  equals  zero  for  all 
values  of  i  and  j  except  for  i=j  since  b5  and  bj  are  indepen¬ 
dent  random  processes,  except  for  the  case  i= j .  Then,  the 
autocorrelation  becomes: 

^(K*TyyT)=<n(*,y)n(xT,yT)>=52  <bi(x)bi<xT)>bd<y-yj)bd(yT-yj)  (A.  47) 

i 

Assuming  all  the  video  lines  are  statistically  equivalent, 
<bi(x>  *s  indePendent  of  i#  and  therefore: 


72 


I 


R(xxTyyT)-<b<x>b(xr>>52  ^1d(y-yj)^1d(yT-y1)  •  (A.  48) 

i 

After  approximating  (A. 48)  by  an  integral  (i.e.,  the  sampling 
artifacts  of  the  display  are  being  ignored) ,  the  resulting 
form  is: 


R<xxTyyT)“<b<x)b(xT>>-2^J^bd<y-yi)bd<yr-yi)dyi  (A.  49) 

where  Ayf  is  the  angular  spread  between  adjacent  vertical 
video  lines. 

This  equation  can  be  recast  by  using  the  appropriate 
change  of  variables: 

p=yT-y5 

Y=y-yT 

X=X-XT 

Assuming  the  process  in  x  is  wide-sense  stationary  (<b(x)b<*T)> 
depends  only  upon  x-xT)  : 


Rm)=R<xxTyyT)“  <b<x)b(xT)>'2^:J-  hd(p>hd(Y<p)dP  *  (A.  50) 

Then,  to  obtain  the  power  spectrum  of  the  noise,  the  Fourier 
transform  of  R(XY)  is  performed: 


'(fx.fy) 


f- 


<b<X)b(XT)>e  X  dX-2y:Hd<fy)H 


d(fy) 


.  Joules2 . 
mrad2 


(A. 51) 


after  use  of  the  Correlation  theorem  [Ref.  39).  The  Fourier 
transform  of  <b(xjb  generates  the  voltage  noise  PSD  (prior 
to  the  display)  "  provided  the  units  are  properly  transformed 
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from  'voltage'  and  'voltage'  space  to  radiant  energy  and 
display  space  [Ref.  30]".  In  other  words,  s(fx,fy)  is  the 
displayed  noise  PSD  and  it  is  expressed  in  radiant  energy 
units/mrad.  Assuming  S(fxify)*S{fjt>S(f  )#  where: 


’(f.) 


-f- 


<b<x)t<xT)>e  dx 


Joules* 

mrad 


(A. 52. a) 


»(fy) 


1  H 

w, a!t 


y'  d(fy>[ 


mrad 


]. 


(A.  52  .b) 


Noting  that  there  is  a  linear  relation  between  voltage  and 
radiant  energy,  S(fx)  can  be  reinterpreted  as  an  electrical, 
optical  or  thermal  noise  PSD  provided  the  appropriate  unit 
conversion  is  performed  to  leave  S(fx  fy)  unchanged.  Applying 
(A. 6)  : 


S0(fx)=Si(fx)HElECT  (fx> 


(A. 53. a) 


f 


or 

s  * 

S0(fx)=r’Q~~“HEiECT  (f.)  tP°wertenp  mrad]  <A.53.b) 

where  r  is  a  constant  such  that  the  expression  can  be  referred 
to  temperature  units  instead  of  volts.  This  is  performed  in 
order  to  make  the  units  of  the  noise  compatible  with  those  of 
the  signal.  In  the  next  paragraph  it  is  shown  that  the 
correct  power  temperature  units  is  ( * C) 2 . 

To  determine  the  value  of  this  constant  r,  the  concept 
from  (A. 19)  can  be  used.  If  both  signal  and  noise  are 
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expressed  in  temperature  units,  by  setting  S/N=l,  the  signal 
(NETD  by  definition)  is  equal  to  the  rms  noise  voltage,  and, 
since  f  (temporal  frequency )=vx  (scanning  velocity)  .  f x  (spatial  frequency  in  x)  , 
(A. 20)  can  be  recast  to: 


NETD=a=  [  f"r  _^h|iki(,  dfjMr  Afx„]5-[ 
j°  5  i(frK)  * 


r  A-f  1 

1  n]  2  [-C]  (A. 54) 


where  Af^  is  the  noise  bandwidth  defined  in  terms  of  spatial 
frequencies.  Using  the  more  commonly  seen  reference  bandwidth 
in  Hz,  the  constant  r  becomes: 


netdV.  , 

r= — __ — U  [  *  c2  mrad ]  (A. 55) 

Afn 

Then,  the  noise  power  spectral  density  (A. 52. a)  in  x  can  be 
represented  by: 


s  ■ 

r  • <b<x>b<Xr)>e12',"*dx  -  NETD2^.v,_i^4ECT  I(j)  f  -c*  mrad]  (».5«) 


where  the  arrow  indicates  an  allowed  substitution  after 
converting  from  voltage  to  thermal  units. 

Since  an  expression  for  the  displayed  noise  power  spectrum 
is  expected,  the  temperature  value  of  the  NETD  must  be 
converted  into  radiant  energy,  following  a  similar  approach  as 
in  (A. 38).  Then: 


NETD'^kliEEP  [ Joules/mrad]  (A. 57) 

v* 

and  after  substitution  of  NETD1  for  NETD  in  (A. 56) 
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Jl<b<x)b(xT>>ej2*fxXdx  - 


k2  NETD2  S><fx)  „2 
"vxAfn  8^*“"  CM 


r  Joules2  n 
L  mrad  J 


(A. 58) 


Thus,  recombining  (A. 51),  and  including  the  display  transfer 
function  along  x  (H^j)  : 


S(fX,fy) 


k2  NETD2  SUM  u2  „2  r  Joules2  t 

^"Afn  <f*>Hd<fx'fy>  [  -^T  3  * 


(A. 59) 


As  previously  stipulated,  expression  (A. 59)  for  the  displayed 
noise  spectral  density  has  the  same  dimensions  as  (A.  51).  From 
(A. 10)  the  matched  filter  can  be  defined  as: 


n«rf(fx.fy)=  (^(fx^Kfyj^tfy))  * 


(A. 60) 


where  Hy^  is  the  normalized  undegraded  Fourier  transform  of 
one  bar  along  x  (See  Appendix  C) .  After  substitution  of 
(A. 59)  and  (A. 60)  into  (A. 37. a),  the  matched- filter  noise 
(perceived  noise)  results: 


Nf-I 


k2NETD2 


vxAy 


;td  r«  r* 

W*  ■*« 0 ' 


r— —  HEIECT  ( f  x  )Hd(  f  x,  f  y  f  x  )HT  ( f  yj^c  f  y)d  f  Xd  f  y 

l(Trx) 


2  (A. 61) 


[  Joule s/mrad2] . 

It  is  noted  that,  as  required  by  the  S/N  calculation  to  be 
discussed  in  the  next  paragraph,  (A. 61)  and  (A. 44)  are  dimen¬ 
sionally  consistent. 

Finally,  the  expression  for  the  MRTD  results  from  combin¬ 
ing  (A. 44)  with  (A. 61)  and  solving  for  AT  with  a  threshold  va¬ 
lue  of  signal-to-noise  ratio  Y.  An  improvement  factor  of  Fr.te 


V 


1 


♦ 
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(frame  rate  x  eye  integration  time)  is  applied  to  the  signal 
since  it  is  assumed  a  process  of  summation  of  the  signal  and 
noise  over  the  frames  in  an  eye  integration  time.  Thus,  the 
MRTD  expression  obtained  is: 

MRTD 


Hoc^x)LfJ.Hmy,^(Vdfv 


,Ayfi&°FrteJlJI S '  j(f  *  jH^ect  cf*)Hd<fx.VH«fx>I^«VH®(,y>dfxdf> 


or 


—  Tnetd 
8 


(A.62) 


v*Ay,-  rr”  UT*> u2  «2  u2  u2  u2  ** 
AfnF7teJ  --JoS'i<f  ,  ELECT  (^x)Hd<fx,fy)HW(fx)HT(fy)M0(fy)arxary 


i<*X>  tt2 


Expression  (A.62)  may  be  somewhat  misleading  since  an  osten¬ 
sible  dependency  of  the  MRTD  on  the  frame  rate  (Fr)  shows  up. 

However,  the  NETD  has  a  similar  dependence  on  the  same 
parameter  (equations  (A. 32)  and  (A. 33))  and  the  overall  effect 
cancels  out.  As  pointed  out  by  Lloyd's  observations 
[Ref.  10,  p. 189] ,  the  MRTD  does  not  depend  on  the  frame  rate. 

It  is  worth  noting  that,  in  the  limit  fox  -*  0  [cy/mrad] , 
the  MRTD  is  actually  lower,  i.e.,  more  optimistic,  than  the 
NETD  (A. 32).  This  is  qualitatively  explained  by  the  intro- 
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duction  of  an  eye-brain  matched  filter.  It  has  been  deter¬ 
mined  that  the  improvement  factor  for  the  sample  system  (App¬ 
endix  D)  is  approximately  1.58.  See  Chapter  III,  Section  B. 


i 
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APPENDIX  B  -  DESCRIPTIVE  FLOWCHART  OP  SIMULATION 


79 


APPENDIX  C  -  EXACT  FOUR-BAR  FOURIER  SPECTRUM 


The  standard  4-bar  pattern  and  the  coordinates  considered 
in  this  work  are  shown  in  the  Figure  C.l  : 


W  -  Angular  width  of  a  bar 

Figure  C.l.  MRTD  Standard  target 

The  bar  pattern  can  be  represented  by  the  following  two- 
dimensional  function: 


it<«.y>  '  t  R«ct 

+  Rect  ]  Rect  (j-\ 


(nr) 


+  Rect 


m 


(C.l) 


The  two-dimensional  Fourier  transform  is  defined  as 


follows: 


£  iT(x.„  e-<  *  Vv>  dxdy  (c.2) 

where  k  =  2irf,  and  k  =  2irf  .  Both  spatial  frequencies,  f  and 

x  a  y  y  a 

fy,  are  expressed  in  cvcles/mrad. 

Since  this  two-dimensional  function  can  be  separated  in 
two  independent  factors  such  that  iT(xy)  =  iT(x)iT(y),  the  two- 
dimensional  Fourier  transform  can  be  expressed  as: 

*20  {  ^T(x,y)  ^  =  *^  ^T(x)  ^  *\0  ^  ^T(y)  ^  (C.3) 

Then: 


*20 


^(x.y)  *  - 


^j0{  Rect 


-j  <kx.x)  ,  2  •  ,,  W. 

e  x  dx  =  sin(kx_) 
k  *2 


(C.4) 


Applying  Shifting  theorem  for  the  four  bars  of  the  target: 


*iD{iT{x)}  =  ^-sin(kx  |)  [e'jkxW  +  eik*w  +  e'j3kxM  +  ej3kxW]  (C.5) 


i.[-  sin(kx|)  +  sin(3kx|)  - 
sin(5kx|)  +  sin(7kx|)  ] 


W[-Sa(kx|)  +  3Sa(3kx|)  - 

5Sa(5kx|)  +  7Sa(7kx|)  ] 


(C.6) 


(C.  7  ) 


where  Sa(u)=sin(u)/(u) .  Similarly,  in  the  y-dimension: 
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(C.8) 


*io  ^  ^T(y)  ^ 


■; 


7“ 

2 

.7“ 

2 


4-j(ky.y) 


dy  = 


2  8111(71^1) 


w. 


=  7  W  Sa(7ky.=  ) 


Finally,  the  two-dimensional  Fourier  transform  of  the 
standard  4-bar  pattern  results: 


•Syitu.y)5  -  7W2[-Sa(kx|)+3Sa(3kx|)- 

5Sa(5k  x|)+7Sa(7kx|)]  Sa(7k") 


(C.9) 


or  expressed  as  a  function  of  f  ,  f  : 

*  7 


IT(fx,fy,  =  ^iT<x,y)>  =  7W2[-Sa(jrfxW)  +3Sa(3?rfxW)  - 

5Sa(5irfxW)+7Sa(7irfxW)  ]  Sa(7rrfyW) 


(C.10) 


which  is  coincident  with  the  expression  employed  by  J.G. 
Vortman  and  A.  Bar-Lev  in  the  model  discussed  in  Chapter  II 
[Ref.  8]. 

The  target  spectrum  is  normalized  to  obtain  the  unitless 
HT(fy)»  HT(f„>'  anc*  nwcfx>  employed  in  the  Ratches  [Ref.  5]  and  Vor¬ 
tman  and  Bar-Lev  models  [Ref.  8].  Thus: 


Ht<V  “  ^nr  =Sa<7kr”>' 


(C.ll) 


w, 


1  W 

HT(f  )=_Li:=  it-Safk^)  +  SSapk,:!)  - 
T(t*J  4  W  4  ^2  ^2  (C.  12) 

5Sa(5kx^)  +  7Sa(7k  ”)],  and 

4u 
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(C. 13) 


}=Sa(kx|)  . 


A  plot  of  the  normalized  two-dimensional  spectrum,  is 
shown  in  Figure  C.2,  where  an  angular  width  of  the  bars  W=1 
mrad  and  a  normalization  factor  of  1/(28W2)  are  used. 

Figures  C.3  and  C.4  show  the  normalized  spectrum  expressed 
as  a  function  of  f1  =f  J  (along  the  x-axis)  and  f'  =f  W  (along 
the  y-axis)  respectively,  for  the  same  angular  width  and 
normalization  factor  as  Figure  C.2.  This  allows  an  estimation 
of  the  two-dimensional  spectrum  for  different  standard  pattern 
sizes. 


Figure  C.2.  Two-dimensional  normalized  target  spectrum 


-1  -D.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 

fx 

Figure  C.3.  Normalized  target  spectrum  along  the  x-axis 
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NORMALIZED  TARGET  SPECTRUM  ALONG  THE  Y  AXIS  (fx=0) 


Figure  C.4.  Normalized  target  spectrum  along  the  y 


APPENDIX  D  -  SAMPLE  SYSTEM 

In  order  to  simulate  the  different  models  of  MRTD  present¬ 
ed  in  this  thesis,  a  sample  system  was  selected  in  coincidence 
with  the  one  employed  by  Lloyd  [Ref.  10].  The  parameters  of 
that  system  are  as  follows: 

i  Lens  focal  length  (d)  = . . . 50000  /xm 

Diameter  of  lens  (D)  = . 20000  /xm 

ft 

F/  number  ( F=d/D)  = . 2.5 

Detector  array  individual  element  size  (square)  = . 0.005  cm 

Horizontal  detector  angular  subtense  (Ax)  = . 1  mrad 

Vertical  detector  angular  subtense  (Ay)  = . 1  mrad 

Detectors  cold  shielding  scheme^ . not  background  limited 

Characteristic  wavelength  of  the  detectors  (lp)== . 11.5  /xm 

Spectral  bandpass  of  detectors^ . 8  /xm  -  11.5  /xm 

Specific  detectivity  at  Xp  (D*Up))  = . 2*1010  cm  Hz1/2/watt 

Specific  detectivity  in  the  bandpass= . lD*(Ap)/lp 

Frame  rate  (Fr)  = . 30  Hz 

Scan  rate  format= . 60  fields/sec 

i  Number  of  detectors  in  parallel  (np)  = . 150 

3  Number  of  scan  lines= . 300 

Inter  lace* . 2  to  1 

Horizontal  scan  efficiency  (r?h)  = . 0.8 

Vertical  scan  efficiency  (rjv)= .  . 0.8 

Overall  scan  efficiency  (n8C=J?h*nv)  = . 0.64 
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Overscan  ratio  ( r]msc)  = . 1 

Distance  between  horizontal  scan  lines  (Ay()  = . 1  mrad 

Horizontal  field  of  view  (HFOV)= . 400  mrad 

Vertical  field  of  view  (VF0V)= . 300  mrad 

Detector  dwell  time  (rd=n*r?sc*X*Y/ (HF0V*VF0V*Fr)  =  .  .2.67*10*5  sec 

Horizontal  scanning  velocity  (vx=X/Td)  = . 37453.2  mrad/sec 

3-dB  frequency  electronic  roll-off  (fr  in  Hertz) =... 18716. 6  Hz 
3  dB  frequency  electronic  roll-off  (frx  in  cy/mrad)=0.5  cy/mrad 


CRT  spot  size  parameter  (a)  = . 1.234 

Noise  equivalent  reference  bandwidth  (Afn)  = . 29.4  Khz 

S/N  threshold  for  detection  of  one  bar  (T)  = . 4.5 

Background  temperature  (TB)  = . 300*K 

Monochromatic  wavelength  of  the  target  ( A. )  — . 10  /xm 

Optical  efficiency  of  the  viewer  (VoM)~ . 0.8 
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